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PREFACE 


This  report  summarizes  the  results  of  the  comprehensive  program  that  has 
been  pursued  to  delineate  all  aspects  of  conductive  interference  in  rail 
transit  systems,  and  to  help  avoid  its  effects  in  new  rail  transit  systems. 
The  program  has  been  undertaken  as  a cooperative  venture  of  U.S.  manufacturers 
of  rail  transit  propulsion  and  signaling  equipment,  rail  transit  system 
operators,  and  members  of  the  research  and  consulting  community,  under  the 
aegis  of  the  Rail  Transit  EMI/EMC  Working  Group,  sponsored  by  the  Urban  Mass 
Transportation  Administration  of  the  U.S.  Department  of  Transportation. 

This  effort  has  enjoyed  the  direct  and  indirect  efforts  of  many  people. 
In  addition  to  Messrs.  Krempasky,  Clark,  Frasco,  Hoelscher,  Rudich,  and 
Truman,  who  are  noted  in  the  References  section  of  this  report,  the  author 
wishes  to  express  gratitude  to  Lennart  Long  of  DOT/TSC,  whose  encouragement, 
support,  and  personal  involvement  were  directly  instrumental  in  the  completion 
of  this  project.  Charles  Edelson,  then  of  SDC,  Inc.,  was  an  active 

participant  in  many  of  the  tasks  reported  on  here.  Donald  Stark  of  Union 
Switch  & Signal  Div.,  American  Standard,  Inc.,  suggested  to  the  author  the 
form  of  the  lumped  equivalent  circuit  shown  in  Sec.  6,  and  provided  much 
useful  information  on  the  nature  of  conductive  interference  in  the  real  world. 
Tatiana  Vinnikova,  then  of  SDC,  Inc.,  was  instrumental  in  the  computation  of 
probability  integrals  in  Section  8.  The  personnel  of  the  Massachusetts  Bay 
Transit  Authority  (MBTA)  gave  invaluable  assistance  for  the  rail  impedance 
tests  described  in  Sections  3-5,  as  did  Irving  Golini,  Michael  West,  John 
Lewis,  and  Paul  Poirier  of  DOT/TSC,  and  Klaus  Frielinghaus  of  GRS,  Inc.  Bay 
Area  Rapid  Transit  District  (BART)  personnel  expended  immense  effort  and 
hospitality  to  make  the  multi -car  measurements  outlined  in  Section  9 possible. 
Final  thanks  go  to  the  sponsor  of  the  UMTA  program  in  rail  transit  EMI/EMC  - 
Ronald  Kangas,  Chief  of  the  Design  Division,  UMTA  Office  of  Systems 
Engineering. 
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EXECUTIVE  SUMMARY 


The  purpose  of  this  report  is  to  summarize  the  results  of  the  comprehen- 
sive program  that  has  been  pursued  to  delineate  all  aspects  of  conductive 
interference  (Cl)  in  rail  transit  systems,  and  to  help  avoid  its  effects  in 
new  rail  transit  systems.  This  program  has  been  undertaken  as  a cooperative 
venture  of  U.S.  manufacturers  of  rail  transit  propulsion  and  signaling 
equipment,  rail  transit  system  operators,  and  members  of  the  research  and 

consulting  community,  under  the  aegis  of  the  Rail  Transit  EMI/EMC  Working 

Group,  sponsored  by  the  Urban  Mass  Transportation  Administration  (UMTA)  of  the 
U.S.  Department  of  Transportation  (DOT). 

In  rail  transit  applications.  Cl  is  defined  as  electromagnetic 
interference  (EMI)  affecting  signaling  or  automatic  train  control  (ATC) 
systems,  caused  by  spurious  currents  flowing  in  the  running  rails  ahead  of  or 
behind  the  train.  A prime  cause  of  Cl,  and  the  one  dealt  with  in  this  report, 

is  the  induction  of  unwanted  currents  in  the  running  rails  by  ac  harmonic 

currents  present  in  the  third  rail.  These  EMI  currents  can  cause  interfering 
voltages  in  wayside  track  circuits  and  carborne  ATC  receivers.  The  ac  harmonic 
currents  in  the  third  rail  can  be  produced  by  rectification  substations  and  by 
solid-state  propulsion  control  systems  on  board  rail  transit  vehicles.  (See 
Ref.  15  for  an  overview  of  rapid  transit  propulsion  systems,  including 
rectification  substations.)  The  major  part  of  the  Cl  program  has  dealt  with 
Cl  arising  from  these  sources. 

Under  the  Cl  program,  results  have  been  obtained  in  the  following  areas: 

• Dependence  of  ac  rail  inductance  on  frequency  and  dc  current. 

• Track  impedances  and  inductive  coupling  between  third  rail  and  track, 
t Analysis  of  the  circuit  consisting  of  third  rail,  track,  track 

ballast,  signaling  components,  and  cars. 

• Statistical  analysis  of  multi -car  Cl. 

t Development  of  Recommended  Practices  outlining  tested  and 

standardized  techniques  for  gathering  Cl  data  and  testing 

compatibility  of  propulsion  and  signaling  equipment. 
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This  report  documents  the  efforts  to  achieve  the  understanding  of  the 
nature  and  characteristics  of  Cl  that  is  required  for  intelligent  application 
of  the  Cl  Recommended  Practices  to  specific  situations. 

As  a result  of  this  program,  rail  transit  Cl  due  to  car  propulsion  systems 
and  rectifier  substations  interfering  with  track  signaling  and  ATC  systems,  is 
well  understood,  predictable,  and  measureable.  The  techniques  of  data 
gathering,  measurement,  and  analysis  developed  in  this  program  and  presented 
in  this  report  should  serve  as  a basis  for  mitigating  the  effects  of  Cl  in 
future  rail  transit  systems. 


1.0  INTRODUCTION 


Rapid  transit  systems  in  the  U.S.  traditionally  have  used  dc  propulsion 
control  systems  relying  on  switched  field  resistors  and  windings.  Their 
signaling  systems  have  been  based  on  power-frequency  track  circuits  and 
wayside  signals.  In  recent  years  the  desire  for  greater  operational 
efficiency  has  led  to  the  introduction  of  advanced  audio-frequency  automatic 
train  control  (ATC)  and  signaling  systems,  and  new  types  of  propulsion  control 
systems.  Advanced  electronics  are  used  in  both  the  new  ATC  systems  and 
propulsion  control  systems. 

Audio-frequency  ATC  systems  can  be  installed  on  continuously  welded  rail, 
thus  obviating  the  need  for  insulated  joints.  In  addition,  such  systems  allow 
easy  integration  of  train  detection  and  transmission  of  ATC  signals  to  the 
train. 

Propulsion  control  systems  using  switched  thyristors  eliminate  the  moving 
parts  and  electrical  contactors  that  prove  to  be  high-maintenance  items  in 
older  dc  propulsion  control  systems.  In  the  future,  the  U.S.  may  see  a 
further  advancement  now  used  in  Europe  and  Japan  - ac  induction  motors  with 
electronically  controlled  dc-to-ac  converters.  These  will  replace  the  dc 
traction  motors  and  their  commutator  rings  and  brushes.  [Ref.  15] 

To  realize  the  advantages  these  modern  electronic  systems,  their  electro- 
magnetic compatibility  (EMC)  must  be  assured.  Early  U.S.  tests  showed  that 
propulsion  systems  with  solid-state  control  can  cause  audio-frequency  harmonic 
currents  to  flow  in  the  third  and  running  rails  and  into  track  circuit 
receivers.  The  problems  of  this  electromagnetic  interference  (EMI)  can  be 
solved,  but  solution  requires  thorough  understanding  of  the  mechanisms 
involved,  and  a coordinated  design  of  the  ATC  and  propulsion  systems. 

The  U.S.  Department  of  Transportation's  EMI/EMC  program  has  the  objective 
of  assuring  the  EMC  of  rail  transit  electrical  and  electronic  subsystems. 
Major  goals  of  this  program  are  to: 
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• Develop  and  validate  Recommended  Practices  (RP's)  for  measuring  EMI  and 
susceptibility  of  rail  transit  electrical  and  electronic  subsystems. 

0 Provide  EMI/EMC  support  to  the  UMTA/STARS  New  Technology  Assessment  and 
Development  Program. 

0 Develop  engineering  guidelines  for  transit  authorities  and  suppliers, 
to  aid  in  the  evaluation  and  elimination  of  EMI  that  could  cause  safety 
or  operational  problems. 

0 Coordinate  and  support  the  activities  of  the  International  Rail  Transit 
EMI/EMC  Technical  Working  Group  (TWG),  and  support  the  APTA/EMI  Liaison 
Board. 

The  purpose  of  the  RP's  is  to  have  a set  of  test  procedures  with  industry- 
wide acceptance  that  can  provide  a common  framework  for  eliminating  and 
avoiding  EMI  problems  and  assuring  EMC.  To  date,  two  sets  of  RP's  have  been 
developed  and  have  been  released  to  the  IEEE  for  their  consideration  as  new 
standards  [Ref.  12].  These  sets  deal  with  two  specific  types  of  EMI  - 
Inductive  and  Conductive. 

As  defined  in  rail  transit,  inductive  EMI  results  from  stray  ac  magnetic 
fields  from  a car's  propulsion  or  auxiliary  power  system  inductively  coupling 
to  the  conducting  loops  under  the  car  comprised  of  signal  receiver  leads, 
rails,  and  car  axles.  Conductive  EMI,  the  topic  of  this  report,  results  from 
the  injection  of  ac  harmonic  currents  from  substation  or  car  into  the 
conducting  loop  formed  by  the  third  rail  and  running  rails,  and  thence  into 
signal  receiver  leads  attached  to  the  running  rails. 
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2.  THE  CONDUCTIVE  INTERFERENCE  MODEL 


This  section  summarizes  the  analysis  of  conductive  interference  in  rapid 
transit  systems.  Figure  2.1  shows  the  configuration  of  third  rail,  running 
rails,  train,  and  substation  that  can  lead  to  the  induction  of  Cl  currents 
into  the  running  rails.  This  figure  shows  balanced  2-rail  audio-frequency 
track  circuits  with  continuously  welded  rail.  The  third  rail  is  positioned 
assymetrically  with  respect  to  the  running  rails,  leading  to  non-zero  mutual 
inductance  between  the  third  rail  and  the  loop  formed  by  the  running  rails. 
If  an  overhead  catenary  were  used  to  carry  dc  current  to  the  train,  and  if  it 
were  positioned  exactly  above  the  centerline  of  the  running  rails,  the  mutual 
inductance  from  it  to  the  running  rail  loop  would  be  zero. 

Cl  currents  induced  into  the  running  rails  do  not  necessarily  result  in 
currents  into  track  circuit  receivers  shunting  the  running  rails.  In  the 
example  in  Figure  2.1,  if  the  third  rail  runs  continuously  down  one  side  of 
the  track,  equal  but  opposite  currents  through  a track  circuit  receiver  will 
flow  due  to  currents  induced  into  the  loops  to  the  left  and  to  the  right  of 
the  receiver,  giving  a net  Cl  current  of  zero.  However,  if  the  third  rail  is 
switched  from  one  side  of  the  track  to  the  other  at  the  location  of  a track 
circuit  receiver,  the  Cl  current  in  the  track  circuit  receiver  can  be  as  great 
as  twice  the  circulating  or  di f ferential -mode  Cl  current  in  the  rails. 
Therefore,  the  layout  of  the  third-rail  circuit  can  have  a direct  impact  on  Cl 
—a  fact  that  track  designers  should  be  aware  of. 

A detailed  analysis  of  this  circuit  must  also  include  the  impedances  of 
track  circuit  components,  and  the  leakage  conductance  of  the  track  ballast. 
As  will  be  shown  in  Section  6,  a good  rule-of-thumb  for  the  ratio  of 
third-rail  ac  current  to  induced  differential  mode  current  in  the  running 
rails  is  the  ratio  of  self-inductance  per  meter  of  the  running  rail  loop  to 
the  mutual  inductance  per  meter  between  third  rail  and  the  running  rail  loop. 
This  ratio  is  approximately  6:1. 

Substations  serve  as  one  source  of  ac  current  in  the  third  rail. 
Substations  generally  behave  as  practically  ideal  voltage  sources  with  very 
low  source  impedance.  Three-phase  and  six-phase  full -wave  rectifier 
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substations  generally  produce  their  strongest  harmonics  at  multiples  of  6 and 
12  times  the  power  frequency  respectively.  However,  all  harmonics  of  the 
power  frequency  can  be  present  in  some  amount  - sometimes  a large  amount  when 
assymetries  exist  in  a substation.  Harmonic  voltage  and  current  content 
varies  as  the  dc  current  drawn  from  the  substation  varies. 

Transit  vehicles  that  employ  chopper  propulsion  control  produce  ac  current 
with  components  at  harmonics  of  the  chopper  frequency.  Generally,  the  chopper 
frequency  is  well-defined.  However,  some  chopper  propulsion  control  systems 
employ  frequency-sweeping  or  pulse-skipping  for  starting  vehicles  from  rest  or 
smoothing  the  transition  to  maximum  power,  resulting  in  production  of  current 
at  frequencies  between  the  normal  harmonics. 

As  will  be  noted  in  Section  7,  a multi -car  train  composed  of  N chopper 
cars  can  be  represented  approximately  as  a parallel  combination  of  independent 
current  sources  in  parallel  with  a source  inductance  whose  value  is  the  line 
filter  inductance  of  one  car  divided  by  N.  Of  special  interest  is  the  fact 
that  at  any  single  chopper  harmonic  frequency,  the  contributions  of  the  N cars 
add  randomly  in  phase.  Since  the  chopper  frequencies  of  separate  cars  in  a 
train  are  very  nearly  but  not  exactly  equal,  their  relative  phases  will  be 
approximately  constant  for  the  short  time  required  for  a train  to  accelerate 
from  rest  to  maximum  speed.  However,  in  the  longer  time  intervals  between 
starts,  the  relative  phases  will  have  time  to  re-randomize. 

The  overall  circuit  in  Figure  2.1  appears  as  a conductor  comprised  of  the 
third  rail  in  parallel  with  another  conductor  comprised  of  the  running  rails, 
with  a train  serving  as  a signal  source  at  one  end,  and  a substation  serving 
as  a signal  source  at  the  other  end.  The  magnitude  of  ac  harmonic  current 
depends  strongly  on  the  total  series  loop  inductance  of  this  conducting  pair. 
The  loop  inductance  per  meter  of  the  third-rail,  running-rail  circuit  is 
approximately  0.9  microhenries.  Since  typical  chopper  cars  have  line  filter 
inductances  of  0.5  to  1 millihenry,  the  inductance  of  the  third-rail,  running 
rail  circuit  is  actually  larger  than  the  net  inductance  of  an  N-car  train,  for 
a train-to-substation  distance  greater  than  500/N  to  1,000/N  meters. 


-5- 


3.  DEPENDENCE  OF  AC  RAIL  IMPEDANCE  ON  FREQUENCY  AND  DC  CURRENT 


3.1  Electromagnetic  Properties  of  Steel  Rails 

To  calculate  expected  values  of  harmonic  currents  and  voltages  in  the 
circuit  of  Figure  2.1,  the  impedance  characteristics  of  the  circuit  must  be 
known.  These  characteristics  depend  upon  the  self  and  mutual  inductances  of 
the  rails  due  to  ac  magnetic  flux  outside  the  rails,  and  upon  the  contribution 
to  rail  impedance  due  to  rail  resistance  and  due  to  ac  magnetic  flux  beneath 
the  surface  of  the  rails.  Since  the  rails  are  ferromagnetic,  rail  inductance 
and  resistance  are  functions  of  dc  current  as  well  as  frequency. 

The  ac  magnetic  field  beneath  the  surface  of  the  rail  decreases 
approximately  exponentially  with  depth  into  the  rail,  with  an  attenuation 
coefficient  defined  as  the  skin  depth.  The  magnetic  permeability  of  steel  is 
many  times  greater  than  that  of  aluminum  or  copper.  However,  steel's 
conductivity  is  many  times  less,  and  therefore  the  skin  depth  in  steel  and 
aluminum  or  copper  are  roughly  comparable.  However,  the  presence  of  ac 
magnetic  fields  within  the  high-permeability  steel  leads  to  an  increase  in 
total  circuit  inductance.  The  contribution  to  inductance  due  to  ac  magnetic 
fields  within  the  steel  will  be  referred  to  here  as  the  "internal  inductance." 

Since  skin  depth  is  smaller  when  magnetic  permeability  is  greater,  the 
audio- frequency  resistance  of  steel  rail  is  significantly  greater  than  that  of 
aluminum  or  copper  conductors  of  equivalent  crossection. 

3.2  Measurement  of  AC  Rail  Resistance  and  Inductance 


A series  of  measurements  was  conducted  to  determine  the  variation  of 
internal  small-signal  ac  inductance  of  railroad  rail  as  a function  of 
frequency  and  dc  current.  The  frequency  range  was  from  25  Hz  to  3 kHz,  and 
the  dc  current  range  was  from  0 to  approximately  1,000  amperes.  This  range  of 
dc  current  is  broader  than  that  previously  investigated  by  Holland-Moritz,  who 
made  similar  measurements  [Ref.  1],  The  configuration  of  the  circuit  employed 
is  shown  in  Fig.  3.1.  This  circuit  is  an  adaptation  of  one  used  by  Kennel ly, 
Achard  & Dana,  to  investigate  variation  of  large-signal  ac  inductance  at  power 
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frequencies  [Ref.  2].  The  difference  lies  in  the  addition  of  the  upper  loop 
in  the  voltage  sensing  conductor  pair,  resulting  in  a nulling  circuit. 

In  the  circuit  shown  in  Fig.  3.1,  the  dc  current  source  was  a parallel 
bank  of  storage  batteries  from  motor  buses.  DC  current  was  controlled  by  a 
bank  of  switched  brake  resistors  from  rail  transit  cars.  The  ac  coupling 
transformer  secondary  winding  was  a solenoidal  motor  smoothing  reactor  from  a 
rapid  transit  car.  The  primary  winding  was  formed  by  wrapping  20  turns  of  No. 
20  AWG  wire  around  the  motor  smoothing  reactor. 

A servo-type  Hall  Effect  electronic  current  transformer  (see  Sec.  9.2)  was 
employed  to  simultaneously  measure  dc  and  ac  current  in  the  rail.  The  turns 
ratio  of  the  current  transformer  was  5,000:1,  and  a 10-ohm  sensing  resistor 
was  used,  resulting  in  a transresistance  ratio  V$ense/Is  = (1/500)  ohm. 

A digital  voltmeter  was  used  to  observe  the  dc  component  of  Vsense.  The 

ac  signal  source  was  a digitally  controlled  sine-wave  oscillator  of  extreme 

frequency  stability.  A two-channel  FFT  network  analyzer  was  used  in  the 

transfer  function  mode  to  observe  the  amplitude  and  phase  of  VQ/Vsense  at 

specific  frequencies.  Measurement  frequencies  were  chosen  that  corresponded 

exactly  to  FFT  frequency  lines.  The  FFT  analyzer  served  as  a two-channel 

digital  filter  for  providing  adequate  signal-to-noise  ratio.  Signal-to-noise 

ratio  for  Y and  Y was  directly  observed  by  using  the  FFT  analyzer  in  the 
o sense 

single-channel  mode  to  compare  the  amplitude  of  VQ  and  Vs  to  amplitudes  of 
interference  present  at  other  displayed  FFT  frequencies.  In  all  cases, 
signal-to-noise  ratio  was  better  than  20  db. 

The  circuit  of  Fig.  3.1  was  used  in  the  following  manner:  First,  a length 

of  copper  pipe  with  1-5/8  inch  outer  diameter,  serving  as  the  reference 
conductor,  was  placed  in  the  position  where  the  rail  was  to  go,  referred  to 
hereafter  as  the  test  conductor  position.  Then,  ac  current  was  applied  to  the 
loop  containing  the  copper  pipe,  and  spacing  d'  between  the  return  conductor 
and  the  nulling  conductor  was  adjusted  to  provide  the  lowest  possible 
magnitude  of  induced  voltage  VQ,  at  a frequency  of  1,000  Hz. 
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FIGURE  3.1.  The  circuit  for  measuring  impedance 
properties  of  rail. 
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The  copper  pipe  had  a wall  thichness  much  less  than  its  radius.  Therefore 
to  a good  approximation,  the  current  in  the  copper  reference  conductor  can  be 
considered  a sheet  current.  An  analysis  of  the  geometry  of  the  sensing  loop 
and  the  magnetic  flux  field  of  the  current  loop  shows  that  the  net  flux 
linkages  of  the  sensing  loop  is  exactly  zero,  provided  that  d'  = r ^ . 

Therefore,  in  the  absence  of  ac  resistance  in  the  copper  pipe,  it  would  have 
been  possible  to  make  VQ  exactly  zero,  when  the  copper  pipe  reference 
conductor  was  in  place. 

It  was  not  neccesary  to  make  the  actual  spacing  between  return  and  nulling 
conductors  exactly  uniform  over  their  length,  since  VQ  only  depends  on  total 
flux  intercepted,  and  not  its  distribution.  However,  it  was  necessary  to 
leave  the  positions  of  return  and  nulling  conductors  undisturbed  after  the 
best-case  null  for  VQ  was  obtained. 

After  the  above  circuit  adjustment,  the  copper  pipe  was  removed  and 

replaced  by  40-foot  lengths  of  railroad  rail  of  a number  of  crossections. 

When  another  test  conductor  replaces  the  reference  conductor,  the  magnetic 

flux  linkages  passing  through  the  sensing  loop  change,  due  to  the  different 

size,  crossecti onal  shape,  and  composition  of  the  new  conductor.  The 

inductive  part  of  VQ  is  due  to  the  difference  between  the  reactances  of  the 

test  conductor  and  reference  conductor.  The  resistive  part  of  V is  due  to 

o 

the  ac  resistance  of  the  test  conductor. 

Consideration  of  the  geometry  of  the  circuit  shows  that  when  the  test 
conductor  replaces  the  reference  conductor,  the  only  place  where  magnetic  flux 
patterns  change  significantly  is  very  near  the  replaceable  conductors,  where 
magnetic  flux  has  cylindrical  symmetry.  Since  points  1 and  2 are  many 
conductor  diameters  in  from  the  ends  of  the  test  or  reference  conductor,  end 
effects  are  avoided.  Due  to  the  cylindrical  symmetry,  we  can  relate  VQ  to  the 
difference  in  conductor  inductance,  and  the  rail  resistance: 

Vo  " ^^meas  + jXmeas^  " ^^meas  + (3.1) 

where  I = ac  source  current,  R = test  conductor  resistance,  and  al  is  the 
s meas 

difference  between  the  internal  inductance  of  the  rail  and  that  of  the 
reference  conductor  for  the  9.14  meter  test  length. 
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It  is  important  to  observe  in  Eqn.  3.1,  that  VQ  depends  on  the  difference 
between  internal  inductances  of  the  rail  and  the  copper  pipe,  and  not  the 
value  of  either  itself.  Furthermore,  the  term  internal  inductance  implies 
inductance  due  to  magnetic  flux  inside  some  reference  surface,  as  opposed  to 
the  external  inductance  due  to  magnetic  flux  outside  the  surface.  We  can 
change  the  values  of  internal  inductance  of  both  the  reference  conductor  and 
test  conductor  by  choosing  a new  reference  surface;  however,  this  does  not 
change  their  differences.  As  a matter  of  convenience,  we  have  defined  the 
internal  inductance  of  the  reference  conductor  to  be  zero. 

As  will  be  seen  in  the  data  to  follow,  a|_  will  be  positive  at  very  low 
frequencies,  due  to  internal  magnetic  fields,  but  negative  at  very  high 
frequencies  where  the  skin  effect  excludes  magnetic  fields  from  the  interior 
of  the  rail.  As  noted  above,  this  difference  is  due  to  magnetic  field 
behavior  within  a few  test  conductor  diameters  of  the  test  and  reference 
conductors.  As  is  discussed  in  Section  3.3  below,  the  value  of  this 
difference  can  be  used  in  conjunction  with  the  known  properties  of  copper  pipe 
to  directly  calculate  the  inductance  per  meter  of  railroad  track  composed  of 
two  running  rails. 

3.3  Rail  Resistance  and  Inductance  Data,  and  Track  Impedance 

Using  the  circuit  of  Fig.  3.1,  data  were  collected  for  85  lb/yd  running 
rail,  100  Ib/yd  running  rail,  and  150  lb/yd  NYCTA  third  rail.  Data  for  dc 
currents  between  zero  and  appx.  1,000  amperes  and  frequencies  between  25  and 
3160  Hz  are  shown  in  Tables  3.1,  3.2,  and  3.3.  Included  in  Tables  3.1  and  3.2 
are  values  of  total  track  resistance  and  inductance  per  meter,  calculated  on 
the  basis  of  the  the  measured  impedance  properties  of  individual  pieces  of 
rail . 

The  relations  between  the  parameters  listed  in  Tables  3. 1-3. 3,  that  are 
derived  from  the  transfer  function  H = (v0/vSense^’  are  as  fol^ows: 


Z. 


me  as 


(3.2) 


X 


meas 


= ojAL  = 2irf  AL 
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w = 9.14  meters  = 30.0' 


AL  " \neas/27Tf  ^Lint,rail  " Lint,cu^w 

(3.2  cont'd) 

Rrail  * (Rtneas/9-14)  °hms/mete'- 

R..  ~ 2Rrail  = 0.219Rmaac  ohms/meter 
tk  meas 

A railroad  track  is  formed  by  two  parallel  pieces  of  rail,  forming  a 
structure  that  looks  like  a two-conductor  transmission  line.  The  transmission 
line  series  resistance  per  meter  of  the  track,  R$er,  is  twice  the  resistance 
per  meter  of  a single  rail.  The  transmission  line  series  inductance  per 
meter,  L$er,  can  be  determined  from  the  data  for  AL  as  described  below: 


Figure  3.2  shows  a section  of  track  and  an  equivalent  transmission  line. 
L k 1S  comprised  of  a part  due  to  magnetic  flux  passing  through  the  air 
between  the  rails,  and  a part  due  to  magnetic  flux  internal  to  the  rails.  In 
other  words,  Ltk  depends  partly  on  the  geometry,  and  partly  on  the  composition 
of  the  conductors.  Where  L ^ cu  = Ltk  for  a simulated  track  constructed  with 
thin-walled  copper  pipes  spaced  the  same  as  the  steel  rails,  the  following 
relation  holds:  [Ref.  3] 


"tk,cu 


• <“0/'')cosh'1(d12/2rcu) 


= ( uQ/ir)  1 °ge( d12/rcu)  = 1-617  nh/meter 


(3.3) 


where  uQ  = magnetic  permeability  of  free  space  = 4*  x 10~^  henries/meter, 
d = distance  between  centerlines  = 1.502  meters  (4*-ll") , and  r = radius 

lu  V*  U 

of  the  copper  pipe  = 0.0264  meters  (0.8125").  (The  log  function  is  used  here 
as  a good  approximation  to  the  cosh  function  for  large  arguments.) 


To  find  Ltk  for  actual  railroad  track,  one  adds  to  the  expression  above 
an  amount  equal  to  twice  the  value  of  al — to  account  for  the  internal 
inductance  of  the  two  conductors--resul ting  in  the  relation: 
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FIGURE  3.2.  The  transmission-line  model  of 

railroad  track. 
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(3.4) 


Ltk  = 1.617xl0'6  + 2AL/w 

= 1.617x10"^  + 0.219AL  henries/meter 

It  is  apparent  from  the  data  in  Tables  3.1  and  3.2  that  although  the 
effect  of  dc  current  on  running  rail  resistance  and  inductance  is  measurable, 
the  overall  effect  on  series  transmission  line  track  inductance  is 
negligible,  at  least  for  dc  currents  less  than  1,000  amperes.  This  is  because 
track  transmission  line  series  impedance  Z$er  is  largely  due  to  the  space 
between  the  rails.  For  the  running  rails  the  effect  of  dc  current  on  rail 
resistance  was  barely  perceptible  at  a frequency  of  55  Hz,  and  therefore  data 
were  not  taken  as  a function  of  current  for  frequencies  above  55  Hz. 

The  data  in  Table  3.3  shows  that  the  effect  of  dc  current  on  the  impedance 
properties  of  third  rail  are  very  pronounced,  even  at  3160  Hz.  At  25  Hz,  the 
ac  resistance  is  nearly  cut  in  half  when  dc  current  increases  from  0 to  1000 
amperes,  and  a sizeable  change  in  inductance  occurs.  This  behavior  is 
discussed  further  in  Section  5. 
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TABLE  3.1 


Rail  Impedance  Data  - 85  lb/yd  Running  Rail 


f 

Hz 

'do 

amps 

Zmeas 

yfl 

Rmeas 

yft 

Y 

meas 

AL 

yh 

Rrail 

yfl/m 

Rtk 

yft/m 

Ltk 

yh/m 

25. 

0. 

670. L 16° 

644. 

185. 

1.18 

70. 

140. 

1.88 

25. 

510. 

680. L 16° 

654. 

187. 

1.19 

71.5 

143. 

1.88 

25. 

775. 

668. L 16° 

642. 

184. 

1.17 

70. 

140. 

1.88 

25. 

840. 

666. L 16° 

640. 

184. 

1.17 

70. 

140. 

1.88 

25. 

920. 

664. L 15° 

641. 

172. 

1.09 

70. 

140. 

1.86 

55. 

0. 

864. L 12° 

845. 

180. 

0.52 

92.5 

185. 

1.73 

55. 

405. 

862. L 12° 

843. 

179. 

0.52 

92. 

184. 

1.73 

55. 

835. 

858. L 11° 

842. 

164. 

0.47 

92. 

184. 

1.72 

65. 

0. 

914. L 10° 

900. 

159. 

0.39 

98.5 

197. 

1.70 

100. 

0. 

1106. £ 4° 

1103. 

77. 

0.12 

120.5 

241. 

1.65 

316. 

0. 

2380. £-32° 

2018. 

-1261. 

-0.64 

221. 

442. 

1.48 

1000. 

0. 

8240.^,-59° 

4244. 

-7063. 

-1.12 

465. 

929. 

1.37 

3160. 

0. 

27800. £-69° 

9963. 

-25954. 

-1.31 

1090. 

2180. 

1.30 
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TABLE  3.2 


Rai 1 Impedance  Data  - 100  lb/yd  Running  Rail 


f 

Hz 

amps 

^meas 

un 

^meas 

uo 

Y 

meas 

yft 

AL 

vh 

Rrail 

uft/m 

Rtk 

uft/m 

Ltk 

uh/m 

25. 

0. 

534. L 21° 

499. 

191. 

1.22 

54.5 

109. 

1.88 

25. 

410. 

534. L 20° 

502. 

183. 

1.16. 

55. 

110. 

1.87 

25. 

680. 

514. z.  20° 

483. 

176. 

1.12 

53. 

106. 

1.86 

25. 

730. 

516. l 20° 

485. 

176. 

1.12 

53. 

106. 

1.86 

55. 

0. 

734. L 11° 

721. 

140. 

0.41 

79. 

158. 

1.71 

55. 

560. 

734. l 11° 

721. 

140. 

0.41 

79. 

158. 

1.71 

65. 

0. 

796. L 9° 

786. 

125. 

0.30 

86. 

172. 

1.68 

100. 

0. 

986. L -4° 

984. 

-69. 

-0.11 

107.5 

215. 

1.59 

316. 

0. 

2520.4-40° 

1930. 

-1620. 

-0.82 

211. 

422. 

1.44 

1000. 

0. 

886O.4-660 

3604. 

-8094. 

-1.29 

395. 

789. 

1.33 

3160. 

0. 

29800.4-74° 

8214. 

-28646. 

-1.44 

899. 

1798. 

1.30 
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TABLE  3.3 


Rai 1 Impedance  Data  - 150  Ib/yd  Third  Rail 


f 

i — i 

Cl 

o 

Zmeas 

^eas 

^meas 

AL 

Rrail 

Hz 

amps 

yfl 

uh 

yft/m 

25. 

0. 

708. L 23° 

651. 

276. 

1.76 

30.2 

25. 

550. 

436. L 18° 

415. 

135. 

0.86 

45.4 

25. 

690. 

394. L 17° 

377. 

115. 

0.73 

41.3 

25. 

875. 

364. L 14° 

353. 

88.1 

0.56 

38.6 

25. 

1000. 

358.4  14° 

347. 

86.6 

0.55 

38.0 

25. 

1085. 

348.4  13° 

339. 

78.3 

0.50 

37.1 

55. 

0. 

934.4  12° 

914. 

194. 

0.56 

100. 

55. 

550. 

6OO.4  2° 

600. 

20.9 

0.06 

65.7 

55. 

680. 

560.4  -1° 

560. 

-9.8 

-0.028 

61.3 

65. 

0. 

992.4  9° 

980. 

155. 

0.38 

107. 

65. 

545. 

66O.4  -2° 

660. 

-23. 

-0.06 

72.2 

65 

675. 

614.4  -5° 

612. 

-54. 

-0.13 

67.0 

65. 

1060. 

564.4-12° 

552. 

-117. 

-0.29 

60.4 

100. 

0. 

1194.4  -2° 

1193. 

-42. 

-0.07 

130. 

100. 

550. 

848.4-16° 

815. 

-234. 

-0.37 

89.2 

100. 

680. 

800.4-16° 

769. 

-221. 

-0.35 

84.1 

100. 

1050. 

774.4-26° 

696. 

-339. 

-0.54 

76.1 

316. 

0. 

2500.4-37° 

1997. 

-1505. 

-0.76 

218. 

316. 

545. 

2340.4-51° 

1473. 

-1819. 

-0.92 

161. 

316. 

680. 

2300.4-53° 

1384. 

-1837. 

-0.93 

151. 

316. 

1045. 

2340.4-59° 

1205. 

-2006. 

-1.01 

132. 

1000. 

0. 

8300.4-63° 

3768. 

-7395. 

-1.18 

412. 

1000. 

545. 

8080.4-71° 

2631. 

-7640. 

-1.22 

288. 

1000. 

675. 

7480.4-67° 

2923. 

-6885. 

-1.10 

319. 

3160. 

0. 

29400.4-76° 

7113. 

-28526. 

-1.44 

778. 

3160. 

545. 

28400.4-78° 

5905. 

-27778. 

-1.40 

646. 

3160. 

680. 

29400.4-78° 

6113. 

-28758. 

-1.45 

669. 
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4.  COMPARISON  OF  CALCULATED  AND  MEASURED  TRACK  IMPEDANCE 


The  track  series  resistance  and  inductance  calculated  in  Section  3 on  the 
basis  of  measured  rail  properties  agree  very  closely  with  results  of 
measurements  of  series  impedance  of  actual  track.  The  track  on  which 
measurements  were  made  had  115  lb/yard  rail,  and  thus  was  slightly  larger  in 
crossectional  area  than  the  100  Ib/yard  rail  used  in  rail  impedance  tests. 
However,  since  perimeter  varies  as  (area)1^,  the  perimeter  of  the  115  Ib/yard 
rail  was  only  an  estimated  7%  larger  than  that  of  the  100  lb/yard  rail. 

Measurements  of  series  track  impedance  of  a length  of  rapid  transit  track 
213  meters  (700* ) long  were  made.  The  track  was  composed  of  115  lb/yd  running 
rail  on  wooden  ties  on  crushed  aggregate  ballast.  The  the  running  rails  were 
shorted  at  one  end,  and  an  audio-frequency  signal  source  was  connected  to  the 
open  end.  An  oscilloscope  was  used  to  measure  relative  amplitudes  and  phases 
of  current  and  voltage  waveforms,  at  a number  of  specific  frequencies,  and  the 
short-circuit  input  impedance  Z was  determined  at  each  of  these  frequencies. 

J V 

Ballast  resistance  was  unknown,  and  was  not  measured.  However,  the 
ballast  was  "clean,"  since  the  track  was  part  of  a little-used  siding,  and  it 
is  estimated  that  ballast  resistance  was  in  excess  of  15,000  ohm-meters  (49.2 
ohm-kft,  i.e.,  49.2  "ohms  in  a thousand  feet").  To  assess  our  ability  to 
infer  series  track  impedance  from  measurement  of  Z , we  briefly  must  review 

J V 

transmission  line  theory  [Ref.  4]. 

Railroad  track  essentially  is  a transmission  line.  A transmission  line  is 
characterized  by  series  impedance  Zser  per  unit  length  and  shunt  admittance 
Y$h  per  unit  length.  These  parameters  determine  the  transmission  line's 
propagation  coefficient  and  characteristic  impedance: 


1 /? 

Y = *ZserYsh*  = ° + J'e  per  meter 


(4.1) 


Current  and  voltage  waves  travel  along  a transmission  line  with  spatial 
variation  given  by 

V+(z,u))  = A exp(-yz) 

I+(z,u>)  = v+(z,w)/Zo 

(4.2) 

V (z,w)  = B exp(+Yz) 

I_(z,U))  = -V_(z,(o)/Zo 

The  "+"  and  subscripts  denote  waves  travelling  in  the  +z  and  -z  directions 
respectively.  The  propagation  coefficient  y = a + je  gives  the  spatial  rates 
of  attenuation  and  phase  shift  of  waves  on  the  transmission  line.  For 
lossless  transmission  lines,  a = 0,  and  y is  purely  imaginary.  In  the  case  of 
railroad  track,  non-zero  rail  resistance  and  non-infinite  ballast  resistance 
both  cause  a to  be  non-zero.  Figure  4.1  depicts  waves  travelling  in  the  +z 
and  -z  directions. 


Treating  the  track  as  a transmission  line,  with  series  impedance  Z . = 

ser  «,  lk 

Rser  tk  + ^Xser  tk  ohms/  meter  of  length,  and  shunt  conductance  Y$h  tk  = 
( 1/Rbai last^  mhos/meter  °f  length,  the  propogation  coefficient  y for  the  track 
is 


Y 


^Zser,tkZRballast^ 


1/2 


(4.3) 


For  shorted  sections  of  track  with  length  d sufficiently  short,  i.e.,  for 
I y I d <<  1,  the  shorted  section  of  track  will  have  an  input  impedance  Z$c 
approximately  equal  to  Zsef  tkd.  However,  the  approximation  starts  breaking 
down  for  values  of  I Yd  I much  larger  than  0.1.  In  other  words,  when  conduction 
from  one  rail  to  the  other  through  the  ballast  becomes  almost  as  easy  as 
conduction  clear  down  one  rail,  across  the  short,  and  back,  you  can  no  longer 
measure  the  impedance  of  the  rail  loop  accurately  from  one  end,  because  you 
will  be  looking  also  at  some  ballast  shunting  it. 
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FIGURE  4.1.  Waves  on  a transmission  line. 
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At  a frequency  of  900  Hz,  a value  of 


'zsc‘ 


= 1.7  ohms  was  measured. 


giving  a first-approximation  value  of  |Zc_r  tl.  | = 1.7/213  = 0.00798  ohms/m,  a 

ber , uv.  i /p  _ a 

first-approximation  estimate  of  |y|  = (0.00798/15,000)  = 7.29x10  /meter, 

and  a first-approximation  estimate  of  | yd | = 0.155. 


The  above  estimated  value  of  |yd|  is  somewhat  larger  than  would  be 
desired,  if  we  are  trying  to  determine  Z ..  from  a measurement  of  Z . If 

SG i 5 lK  SC 

ballast  resistance  in  fact  is  higher,  the  estimate  is  more  accurate.  We  will 
assume  that  for  frequencies  up 
val id. 


to  1 kHz,  our  determination  of  Z is 

ser,tk 


Data  for  R$er  tk  and  L$er  tk  = Xser  tk/2Tf  are  shown  plotted  in  Figures 

4.2  and  4.3,  for  frequencies  of  25,  60,  100,  316,  and  900  Hz.  Also  shown  are 

the  data  from  Tables  3.1  and  3.2  for  85  lb/yd  and  100  lb/yd  running  rail,  and 

data  given  by  Carey  for  130  Ib/yd  running  rail  [Ref.  5].  Note  that  there  is 

general  qualitative  agreement  between  the  behavior  of  L$er  tk  and  Rser  tk  vs. 

frequency  as  inferred  from  measurement  of  rail  impedance,  and  as  inferred  from 

measurement  of  Z . 

sc 

One  interesting  aspect  of  this  data  is  how  little  variation  there  is  in 
the  inductance  data  as  a function  of  rail  weight.  In  Section  5 below,  we  will 
develop  a characterization  of  rail  based  on  an  effective  radius.  Since  radius 
varies  as  (weight)l/2,  and  since  the  air  contribution  to  inductance  varies 
logarithmically  with  radius,  little  variation  is  expected. 

In  conclusion,  we  have  found  that  the  audio-frequency  impedance  of 
railroad  track  can  be  accurately  calculated,  based  on  measured  values  of  rail 
impedance.  In  Section  6 below,  we  will  use  the  measured  values  of  rail 
impedance,  together  with  knowledge  of  rail  spacing,  to  calculate  the 
distributed  circuit  properties  of  a circuit  comprised  of  two  running  rails 
plus  a third  rail. 
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FIGURE  4.2.  Transmission  line  series  resistance 

for  railroad  track. 


-21  - 


FIGURE  4.3.  Transmission  line  series  inductance 
for  railroad  track. 
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5.  MEASURED  RAIL  IMPEDANCE  VS.  SKIN  EFFECT  THEORY  [Ref.  6] 


5.1  Skin  Effect  Theory  - the  Skin  Depth  Parameter 


It  would  be  desirable  to  have  a simple  model  that  accounted  for  rail 
impedance  as  a function  of  frequency.  For  non-ferrous  conductors  of  regular 
geometry,  such  models  do  exist.  They  take  into  account  the  magnetic 
permeability,  dielectric  constant,  and  conductivity  of  the  conducting  medium. 
For  planar  conductors,  solutions  for  surface  impedance  involve  exponential 
functions,  as  is  discussed  below.  For  circular  conductors,  solutions  involve 
Bessel  functions;  but  empirically,  impedance  behavior  becomes  very  simply 
described  at  very  high  frequencies. 

Railroad  rail  does  not  present  such  a simple  case.  It  has  very  irregular 
crossection,  and  it  is  made  from  ferrous  material.  We  believe  that  a complete 
accounting  of  high-frequency  impedance  of  steel  rails  must  take  magnetic 
hysteresis  into  account,  probably  by  ascribing  an  imaginary  part  to  the 
magnetic  permeability.  Such  an  effort  is  beyond  the  scope  of  this  project. 
What  follows  is  a description  of  the  Skin  Effect,  and  its  application  to 
steel  rails  in  a qualitative  manner. 

At  high  frequencies,  electromagnetic  fields  are  excluded  from  the  interior 
of  good  conductors,  and  ac  currents  only  flow  on  the  surface  of  the 
conductors.  This  phenomenon  is  referred  to  as  the  Skin  Effect.  Current 

density  decreases  exponentially  with  depth  below  the  surface  of  a conductor 
with  plane  surface.  The  skin  depth  parameter  <s  is  equal  to  the  inverse  of 
the  real  part  of  the  propagation  coefficient  for  electromagnetic  waves  within 
the  conductor.  The  propagation  coefficient  r is  given  by  the  relation 

Y = a + je  = [(jwpMff  + jaje)]^^  (5.1) 

In  a metallic  conductor  at  frequencies  less  than  1016  Hz,  a >>  me  and  the 
resulting  relations  for  -y,  a,  and  6 are 
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Y = ( 1 + j ) ( irfya) 


a = 6 = Re[y]  = (Trfpa)^^ 


(5.2) 


where  u = magnetic  permeability  and  a = conductivity  of  the  conductor. 

5.2  Skin  Effect  in  Conductors  of  Circular  Crossection 

For  a solid  conductor  of  circular  crossection,  variation  of  current 
density  with  depth  is  given  by  Bessel  functions  of  complex  argument. 
However,  at  high  frequencies,  where  6 becomes  much  less  than  the  radius  of 
the  conductor,  decrease  in  current  density  vs.  depth  once  again  becomes 
exponential.  For  conductors  of  irregular  crossection  such  as  railroad  rails, 
only  numerical  solutions  for  current  distribution  exist.  However,  beginning 
at  a specific  point  on  the  surface  of  the  conductor,  as  frequency  increases 
and  skin  depth  becomes  much  less  than  the  radius  of  curvature  of  the  surface, 
current  once  again  will  decrease  exponentially  as  a function  of  depth. 

As  frequency  increases  and  the  effective  current-carrying  crossection  of  a 
conductor  decreases,  the  series  impedance  per  unit  length  increases.  The 
series  resistance  and  internal  reactance,  which  is  that  part  of  reactance  due 
to  fields  within  the  conductor,  each  increase  with  increasing  frequency. 

One  surprising  aspect  of  the  skin  effect  phenomenon  is  that  in  a non- 
ferrous  solid  circular  conductor  at  high  frequencies,  the  internal  reactance 
becomes  equal  to  the  series  resistance.  This  fact  can  be  understood  by 
examining  the  behavior  of  electromagnetic  fields  beneath  the  surface  of  a 
plane  conductor  at  high  frequencies.  The  characteristic  impedance  for  plane 
waves  within  a conducting  medium  is 


ZQ  = [jwp/(o  + jwe)]1/2 


(5.3) 


In  a metallaic  conductor,  at  frequencies  below  10^  Hz,  the  resulting  form  of 
Zo  is 
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(5.4) 


ZQ  = (1  + j ) ( irf y/a)^^ 

Since  the  real  and  imaginary  parts  of  1Q  are  equal,  sheet  current  just 
beneath  the  surface  of  the  conductor  will  give  rise  to  exactly  equal  resistive 
and  reactive  components  of  electric  field  at  the  surface  of  the  conductor. 
Thus,  the  resistive  and  reactive  components  of  "internal  impedance,"  due  to 
fields  and  currents  within  the  conducting  medium,  are  equal. 

Figure  5.1  shows  behavior  of  resistance,  internal  reactance,  and  internal 

inductance  per  unit  length  for  a non-ferrous  circular  conductor  as  a function 

of  (r  /<5),  where  r is  the  conductor  radius.  Note  that  values  of  resistance 
c c 

and  reactance  are  normalized  to  the  low-frequency  resistance  per  meter 

R0  = l/o*rc2  (5.5) 

Note  that  < r / s ) is  proportional  to  f*X2.  For  (r  /«)  > 2,  very  good 
approximations  for  resistance,  internal  reactance,  and  internal  inductance 
are: 


R/R0  = 0.25  + rc/26 

R = 0.25Rq  + Rorc/26 

= 0.25RQ  + ( u/a)  1/2(  l/2ir1/2rc)f 1/2  (5.6) 

Xint  ' “L1nt  ' Vc/2S 

L1nt  - (.V.)1/2(l/4.3'2rc)f-1'2 

The  above  relations  show  that  as  far  as  causing  a large  value  of  L.  is 
concerned,  a large  value  of  y has  no  more  effect  than  a small  value  of  a. 
This  surprising  result  is  somewhat  at  odds  with  intuition,  but  it  is  a direct 
result  of  the  propagation  characteristics  of  electromagnetic  waves  penetrating 
the  surface  of  the  conductor,  as  discussed  above. 
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FIGURE  5.1.  Resistance  and  internal  reactance  for  solid 

circular  conductors,  normalized  to  dc  resist- 
ance, and  internal  inductance  normalized  to 
dc  internal  inductance,  vs.  ratio  of  radius 
to  skin  depth.  From  Ramo  & Whinnery  [Ref.  6]. 
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5.3  Skin  Effect  in  Railroad  Rail 


It  is  not  possible  to  construct  a model  for  railroad  rails  based  on  equi- 
valent non-ferrous  conductors  of  circular  crossection  that  accurately  accounts 
for  both  internal  rail  inductance  and  resistance.  The  behavior  prescribed  for 
such  a model  in  Eqn.  5.6  specifies  very  rigorously  the  relation  between  R, 
Xint*  and  Lint*  as  a ^unctlon  °f  frequency.  Analysis  of  data  for  the  rails 
tested  shows  that  the  behavior  of  the  impedance  parameters  as  a function  of 
frequency  does  not  conform  to  such  a model.  One  reason  is  that  magnetic 
hysteresis  loss  leads  to  higher  skin  resistance;  inclusion  of  hysteresis 
effects  is  beyond  the  scope  of  this  report.  However,  such  a model  can  account 
for  the  behavior  of  rail  inductance  alone  as  a function  of  frequency.  Such  a 
model  is  useful  at  audio  frequencies,  where  the  preponderant  part  of  overall 
track  impedance  is  the  inductive  part,  and  where  the  preponderant  part  of  the 
inductance  is  due  to  the  free  space  between  the  rails.  The  relation  of 
measured  rail  data  to  equivalent  circular  models  is  presented  here. 

Based  on  nominal  data  provided  by  Frielinghaus  [Ref.  7],  for  conductivity 
and  permeability  of  steel  running  rail,  namely 

c 

a . = a /12.1  = 4.79x10  mhos  per  meter 

rail  cu  K 

(5.7) 

y = 20yrt  = 20x4ttx107  = 2.51xl0“5  henries/meter 

rail  0 


the  nominal  value  of  skin  depth  at  25  Hz  is  1 cm.  This  distance  is  much  less 
than  the  radius  of  a circle  whose  area  equals  the  crossection  of  railroad 
rail.  However,  it  is  approximately  equal  to  the  thickness  of  the  web  or  lower 
flange. 

We  can  determine  a value  of  effective  radius  from  measured  inductance  vs. 
frequency.  From  Eqn.  5.6,  the  expected  form  for  measured  rail  inductance  per 
meter  as  described  in  Section  3 is 


Lmeas  ^ yo^ir^°^e^rcu^reff  ^ 


+ (y  /a  . , ) 1/2(  1/4tt 3/2r  ~)f"1/2 
rail  rail  eff 


(5.8) 
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where  reff  is  the  radius  of  equivalent  circular  conductor,  and  rcu  is  the 
radius  of  the  copper  pipe  that  served  as  a reference. 


:-l/2 


In  Eqn.  5.8,  it  is  seen  that  as  f » and  f"x/‘-  -►  o, 

frequeny-dependent  term  disappears,  leaving  a term  dependent  only  on  r and 

-1/2  ^u 

reff*  A Plot  Lmeas  vs*  ^ for  the  Wyd  running  rail  is  shown  in 

Figure  5.2.  A best-fit  straight  line  drawn  through  the  data  points  yields  an 

intercept  at  f"1/2  = 0 of  -1.65  Mh,  or  1.65/9.14  = -0.1805  Mh/meter,  and  a 

corresponding  effective  radius  of 


reff  = rcuexP(0*  1805^/2")  = 5.09  cm  = 2 inches  (5.9) 

The  slope  of  the  line  is  14.5  uh  Hzly/2,  which  divided  by  the  9.14  meter 
length  yields  a slope  of  1.59  uh  m_1  Hz1^2.  Taking  the  value  of  r ~ from 
Eqn.  5.9  and  the  value  of  this  slope,  and  using  Eqn.  5.8  to  solve  for  an 
effective  value  of  (u/o)  based  on  this  inductance  data,  we  obtain 

(m/o)  = I6*3r2ff(1.59  x 10"6)2  = 3.23  x 10“12  ohm-henries  (5.10) 

This  value  of  (m/°),  together  with  the  effective  radius,  can  be  used  to 
calculate  rail  internal  inductance  and  track  inductance  over  the  range  of 
frequencies  covered  with  very  good  accuracy.  Of  course,  the  accuracy  is 
enhanced  by  the  fact  that  as  far  as  the  inductance  of  two-rail  track  is 
concerned,  the  preponderance  of  the  inductance  is  due  to  the  free  space 
between  the  rails. 

Note  that  we  have  no  way  to  separately  determine  u ...  or  o from  our 
data,  since  all  our  data  was  taken  at  frequencies  at  which  skin  depth  was  many 
times  less  than  reff.  The  value  determined  here  is  in  qualitative  agreement, 
but  significantly  different  from  the  corresponding  value  calculated  from  the 

1 p 

data  given  in  Eqn.  5.7  of  5.24  x 10  ohm-henries. 

The  impedance  properties  of  the  150  lb/yd  third  rail  show  striking 

variation  over  a dc  current  range  from  0 to  1000  amperes.  The  data  from  Table 

3.3  have  been  used  to  construct  plots  of  L ac.  as  a function  of  f^/2  for  the 

roeas 
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FIGURE  5.2.  AL  vs.  for  100  Ib/yd  running  rail. 
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2.0 


FIGURE  5.3.  AL  vs.  f“^2  for  150  Ib/yd  high-conductivity 
third  rail,  for  three  values  of  dc  current. 
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third  rail,  for  dc  currents  of  0,  550  amperes,  and  1050  amperes.  These  plots 
are  shown  in  Figure  5.3.  A family  of  straight  lines  has  been  drawn,  one  for 
each  value  of  1^  , with  a common  vertical -axis  intercept,  and  with  a quali- 
tatively good  fit  to  the  corresponding  data  points. 

The  lines  intersect  the  vertical  axis  at  -1.7  yh,  or  -0.1860  yh/meter  for 

-6  -6 

the  9.14  meter  length,  and  the  lines  have  slopes  of  1.89  x 10  , 1.40  x 10  , 

and  1.22  x 10"®  irnT^Hz^2,  for  the  corresponding  dc  currents  of  0,  550,  and 

1050  amperes. 

The  value  of  -0.1860  yh/m  yields  an  effective  radius  for  the  third  rail  of 

r ff  = rcuexp(0.1860yQ/2*)  = 5.23  cm  (5.11) 

Using  the  three  values  of  slope  above  in  place  of  the  1.59  x 10"6  value  in 
Eqn.  5.10,  and  using  the  new  value  of  we  obtain  the  following  estimates 

of  (y/°)  for  the  three  currents: 

Idc,  amps  y/CT,  ohm-henries  for  150  lb/yd  third  rail 

0 4.85  xlO"1  2 

550  2.66  x 10“12 

1050  2.02  x 10~12 

We  can  assume  that  the  value  of  o remains  the  same,  but  that  the  effective 

small -signal  value  of  y decreases  as  dc  magnetic  field  saturates  the  metal 
near  the  surface  of  the  rail.  The  drop  in  y with  increasing  dc  current  leads 
to  an  increase  in  skin  depth,  and  a corresponding  decrease  in  rail  resistance 
as  seen  in  the  data  of  Table  3.3. 

The  drop  in  the  internal  inductance  of  the  third  rail  seen  here  is  large 
enough  to  be  significant  as  far  as  conductive  interference  is  concerned.  In 
Section  6,  we  will  see  that  the  total  inductance  of  the  third-rail  running- 
rail  loop  is  has  a nominal  value  based  on  air  inductance  of  0.9  yh/meter  of 
circuit  length.  At  25  Hz,  we  see  a decrease  here  of  0.135  yh/meter  in  third 
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rail  internal  inductance,  and  therefore  a corresponding  decrease  in  the 
third-rail  running  rail  loop  inductance,  as  current  increases  from  0 to  1050 
amperes. 

This  variation  in  y with  current  also  is  an  indication  that  the  third  rail 
could  be  a significant  factor  in  producing  harmonics  and  intermodulation 
products  of  interference  signals  present  at  chopper  and  substation  harmonic 
frequencies.  We  do  not  know  how  significant  the  third  rail  is  at  causing  such 
mixing  in  comparison  with  other  nonlinear  elements  in  the  system,  but  we  can 
see  that  the  running  rails  will  be  practically  free  of  such  effects  in 
comparison  to  third  rail  such  as  that  tested  here. 

5.4  Conclusions 


Railroad  running  rail  of  a variety  of  weights  is  seen  to  have  internal 
inductance  properties  equivalent  to  solid  steel  circular  cylinders  of  radius 
5.09  cm  (2"),  with  a ratio  (y/a)  equal  to  3.23  x 10  ohm-henries.  Since  air 
inductance  is  the  dominating  term  in  the  series  transmission  line  impedance  of 
railroad  track,  a model  for  total  track  inductance  based  on  these  properties 
gives  accurate  results,  and  magnetic  hysteresis  effects  can  be  neglected. 

Third  rail  of  an  alloy  with  a high  copper  content  was  found  to  have 
internal  inductance  properties  that  varied  significantly  with  dc  current.  The 
150  Ib/yd  rail  tested  had  inductance  properties  equivalent  to  a solid  cylinder 
of  radius  5.23  cm  (2.06"),  and  an  effective  small-signal  magnetic  permeability 
that  decreased  with  increasing  dc  current. 

As  will  be  seen  in  Sections  6.9  - 6.11,  the  internal  inductance  of  the 
third  rail  directly  enters  the  expressions  for  series  impedance  of  the  third- 
rail  running-rail  loop.  At  25  or  60  Hz,  the  loop  inductance  per  meter  will 
decrease  by  almost  0.1  yh/meter  as  current  through  the  third  rail  tested 
increases  from  0 to  1,000  amperes.  And,  we  cannot  say,  based  on  the  data 
taken,  how  much  the  decrease  would  be  as  the  dc  third-rail  current  increased 
to  10,000  amperes.  However,  since  this  loop  inductance  is  dominated  by  air 
inductance,  this  variation  should  not  be  of  actual  concern,  especially  at 
audio  frequencies  where  variation  is  less. 
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6.  CIRCUIT  ANALYSIS  OF  TRACK  WITH  THIRD  RAIL,  TRACK  CIRCUITS,  AND  BALLAST 


6.1  The  Model 


The  model  that  accounts  for  the  audio-frequency  circuit  characteristics  of 
track  with  third  rail,  track  circuits,  and  ballast  is  shown  in  Figure  6.1.  In 
this  model,  the  track  forms  a transmission  line  with  series  transmission  line 
impedance  Z$er,  and  shunt  transmission  line  admittance  Y$h  = l/Rba-j ^ 

Rballast  1S  in  ohm-meters»  an<1  1S  equal  to  the  ballast  resistance  in  ohm-kft 
multiplied  by  304.7.  Inclusion  of  ballast  conduction  yields  a more 
complicated  model  than  than  the  case  without  ballast  whose  analysis  is  shown 
by  Hoelscher  & Rudich  [Ref.  8]  for  this  track  geometry,  and  by  Lowe  and  Mel  lit 
[Ref.  14]  for  the  "fourth  rail"  geometry  used  in  the  London  Underground. 


Figures  6.1(a) ,(b) ,(c)  show  three  track  circuit  configurations:  Double- 
rail track  circuit  terminated  by  impedance  bonds;  single-rail  track  circuit 
with  signal  rail  adjacent  to  third  rail;  and  single-rail  track  circuit  with  dc 
return  rail  adjacent  to  third  rail. 

The  main  objective  in  analyzing  this  circuit  is  to  calculate  the  current 
and  voltage  delivered  to  the  track  receiver  due  to  audio- frequency  current  in 
the  third  rail.  The  coupling  from  third  rail  to  track  receiver  is 
characterized  by  the  conductive  interference  current  transfer  function 


H 


cond 


V!3 


6.2  Double-Rail  Track  Circuit 


(6.1) 


Figure  6.2  shows  the  definitions  of  currents  and  voltages  in  the 
double-rail  track  circuit.  I3  is  the  audio-frequency  third-rail  current.  I^ 
is  also  the  common-mode  component  of  current  in  the  running  rails.  We  shall 
neglect  leakage  conductance  and  capacitance  from  third  rail  to  running  rails 
or  ground.  I^  is  then  presumed  to  be  constant  in  x.  Travelling-wave  effects 
in  the  propagating  structure  formed  by  the  third  rail  and  running  rails  are 
therefore  neglected,  and  this  model  will  be  in  error  for  overall  circuits 
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VKd) 


I, (x)=I(x)-I3/2 


a)  Balanced  case. 


I, (x)=I(x)-I3 


I R= I ( d ) 


■»  > vfxl  I ( x ) 

t < I 2 ( X ) =- 1 ( x ) V(x)  


x = 0 


x = d 


b)  With  third  rail  adjacent  to  signal  rail. 


I, (x)=I(x) 

' 1 R=I  ( d ) 


x = 0 1 3 x = d 


c)  With  third  rail  adjacent  to  dc  return  rail. 


FIGURE  6.1.  The  coupled  transmission  line  model  of  running 
rails  and  third  rail. 
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I n ( X ) = I ( X ) - I 3/2 


+ 


1 2 ( x ) = - I (x)  - I 3/2 


V ( x ) 


I3  = const. 


FIGURE  6.2.  Currents  and  voltages  in  the  balanced 

double-rail  track  circuit. 
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longer  than  approximately  0.1  free-space  wavelengths  at  any  particular 
frequency.  For  instance,  at  5 kHz,  (O.lc/f)  = 6 km.  This  potential  error  is 
not  considered  important,  since  substations  generally  are  spaced  closer  than 
this  distance. 

I ( x ) is  the  differential  or  balanced  component  of  current  in  the  running 
rails.  Individual  rail  currents  in  the  running  rails  are 

Ix(x)  = I (x)  - I3/2 

(6.2) 

I2(x)  = -I(x)  - I3/2 

Voltage  between  the  running  rails  is  V(x). 

The  track  circuit  may  be  either  occupied  or  unoccupied.  If  occupied,  then 
the  train  axles  nearest  the  receiver  provide  esentially  a short  circuit, 
causing  = 0.  The  impedance  bonds  are  designed  to  have  very  low  common-mode 

impedance  but  very  high  differential -mode  impedance  ZB1,  ZB2.  We  shall  assume 
that  the  common-mode  bond  impedance  is  zero. 

The  mutual  inductance  Mgrr  per  meter  between  third  rail  and  running  rail 
loop  is  determined  by  the  geometric  arrangement  of  the  rails.  Figure  6.3 
shows  an  end  view  of  the  rails.  Third-rail  current  I3  gives  rise  to  azmuthal 
magnetic  field  = IB/2irr  circulating  around  the  third  rail,  where  r is 
radial  distance  from  the  third  rail.  Integrating  from  d,^  to  d^  gives 

total  magnetic  flux  * per  meter  passing  through  the  running  rail  loop.  Then, 

M3rr  = $/I3  = (m0/2it)1  nCd^/d^)  (6.3) 

6.3  Differential  Equations  for  Circuit 

Figure  6.4  shows  an  incremental  portion  of  the  circuit  dx  meters  long. 
Application  of  Faraday's  Law  here  yields  the  differential  relation 

V ( x ) - V(x+dx)  = -dV ( x ) = [I(x)Z$er  - jwM^ydx  (6.4) 
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FIGURE  6.3.  Magnetic  coupling  between  third  rail  and 

running  rails. 
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FIGURE  6. 


I (x)— I3/2 


(zser/2)dx 


At 


x + dx 


V ( x ) 


Flux  from  1^ 

>Yshdx  (*Z  * /2)dx  V(x+dx) 
ser 

'V 


I ( x ) +1 3/2 


I 


a)  Double-rail  balanced  case. 


I(x)-I 


Aa 


• • 


I (x) 


At 


1 3 

b)  With  signal  rail  adjacent  to  third  rail 


I (x) 


At 


• • 


I(x)+I3 


At 


c)  With  dc  return  rail  adjacent  to  third  rail 


. An  incremental  portion  of  the  third  rail  and 
running  rails.  Only  the  definition  of  running 
rail  currents  changes  from  case  to  case. 
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which  translates  into  the  differential  equation 


dV(x)/dx  = -I(x)Z$er  + jwM3rrI3  (6.5) 

The  corresponding  differential  relation  for  current  is 

I ( x ) - I(x+dx)  = -dl(x)  = Y$hV(x)  (6.6) 

which  yields  the  differential  equation 

dl(x)/dx  = -Y$hV(x)  (6.7) 

Differentiating  Eqn.  6.5  once  again  while  remembering  that  dl3/dx  = 0,  and 
substituting  in  the  relation  for  dl/dx  from  Eqn.  6.7,  we  obtain  the 
second-order  differential  equation 

d2V(x)/dx2  = YshZserV(x)  (6.8) 

Where  the  complex  propagation  parameter  is 

* ■ <Yshzser>1/2  '6‘9> 

we  can  write  a general  solution  for  V ( x ) in  terms  of  a weighted  sum  of 
rightward-  and  leftward-travelling  waves: 

V(x)  = Aexp(-yx)  + Bexp(+Yx)  (6.10) 

Differentiating  Eqn.  6.10,  substituting  into  Eqn.  6.5,  and  using  the  relation 
ZQ  = ^Zser/Ysh)^2  gives  the  relation  for  I(x): 

I (x)  = (l/Z0)[Aexp(-Yx)  - Bexp(+yx)]  + U3/Zser) jwM3rr  (6.11) 
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6.4  An  Equivalent  Lumped  Circuit 


At  this  point  we  could  solve  any  particular  circuit  with  specific 
terminations  at  x = 0 and  x = d by  finding  the  values  of  A and  B that  match 
the  boundary  conditions  for  V(0)/I(0)  and  V(d)/I(d).  We  would  then  know  V(x) 
and  I ( x)  at  every  point  in  the  circuit.  This  is  more  information  than  we  need 
in  most  circumstances.  What  we  really  want  to  know  are  the  values  of  current 
and  voltage  into  the  terminations  at  the  ends  of  the  circuit.  Therefore,  we 
will  follow  the  approach  outlined  below  of  modelling  the  transmission-line 
structure  by  a simple  lumped-circuit  two-port  network  model  that  contains  an 
internal  voltage  source  to  account  for  voltages  induced  by  the  third  rail 
current. 


The  transmission  line  structure  above  is  symmetrical  end-for-end.  It 
therefore  can  be  modelled  with  an  end-for-end  symmetrical  n-circuit  as  shown 
in  Figure  6.5.  To  find  the  corresponding  relations  for  and  Z2  in  the 

equivalent  IT-circuit,  we  need  only  ask  that  the  transmission  line  and  the 
n-circuit  both  give  equal  input  impedance  Zi  = V(0)/I(0)  looking  in  at  x = 0 
for  two  cases:  Output  end  open-circuited,  i.e.,  1(d)  = 0;  and  output  end 

short-circuited,  i.e.,  V ( d ) = 0.  To  find  Z^^  and  Z2,  we  do  this  with  Ig  = 0. 
We  will  then  let  Ig  be  non-zero  to  find  the  proper  form  of  Vind,  which 
accounts  for  induced  currents  and  voltages  in  the  circuit  due  to  Ig. 

Simultaneous  solution  of  Eqn's  6.10,  6.11  for  the  short-circuit  case  gives 

zisc  = [V(0)/I(0)]v(d)=0  = Z0sinh(Yd)/cosh(Yd)  = Zj_ | |Z2  (6.12) 

where  "||"  is  the  "parallel"  operator:  Za||Zb  = (Za_1  + Zb-1)-1. 


Simultaneous  solution  of  Eqn's  6.10,  6.11  for  the  open-circuit  case  gives 
zioc  = [V(0)/K0)]1(d)  = 0 

(6.13) 

= ZQcosh(Yd)/sinh(Yd)  = Zx | I ( Zj  + Z2) 
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FIGURE  6.5.  The  lumped-element  pi-circuit  that  is 
equivalent  to  the  transmission  line  formed  by  the 
running  rails. 


FIGURE  6.6.  The  lumped-element  pi-circuit  with 
coupling  from  third  rail  included. 


FIGURE  6.7.  The  lumped-element  pi-circuit  with 
third-rail  coupling  and  track  transmitter  and  receiver 
impedances  added. 
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Then  simultaneous  solution  of  Eqns  6.12,13  gives 


= ZQsi nh(Yd)/[cosh(Yd)  - 1] 

(6.14) 

Z = Z sinh(yd) 

2 o 

The  first  step  in  finding  the  value  of  the  induced  voltage  shown  in 
Figure  6.5  is  to  find  the  Thevenin-equi valent  open-circuit  output  voltage  Vth 
that  appears  at  x = d when  the  track  is  shorted  at  x = 0 causing  V(0)  = 0. 
Figure  6.6  shows  the  lumped  circuit  with  Vth  included.  Setting  V ( 0 ) = 0 in 
Eqn.  6.10  yields  the  relation  B = -A  for  this  case.  Using  this  result  to 
eliminate  B from  Eqn.  6.11,  setting  1(d)  = 0 in  Eqn.  6.11,  and  solving  the 
resulting  relation  for  A yields 

A = -B  = -I3ju)M3rr/2YCOSh(Yd)  (6.15) 

Examination  of  the  circuit  shown  in  Fig.  6.6  and  Eqn.  6.10  shows  that 

vth  = vindzl/(zl  + V = v(d)  (6-16) 

Using  Eqn's  6.10,15,16  to  solve  for  Vj  ^ yields  the  result 

Vind  = ^ u>M3rrI3si  slh(  (6.17) 

6.5  Single-Rail  Track  Circuits 

In  the  single-rail  track  circuit  with  the  signal  rail  adjacent  to  the 
third  rail  shown  in  Fig.  6.1(b),  the  signal  current  is  noted  as  I ( x) . An 
analysis  of  an  incremental  portion  of  this  circuit  yields  the  differential 
rel ati on 

V(x)  - V(x+dx)  = — dV ( x ) = [I(x)Zser  - - j^^Ddx 

(6.18) 

■ Cl(x)Zser  - + Zser/2)]dx 
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Writing  this  as  a differential  equation,  we  have 


dV(x)/dx  = -I(x)Zser  + I3(j“M3rr  + Z$er/2)  (6.19) 

In  the  single-rail  track  circuit  with  the  dc  return  rail  adjacent  to  the 
third  rail  shown  in  Fig.  6.1(c),  the  signal  current  is  also  noted  as  I ( x) . An 
analysis  of  an  incremental  portion  of  this  circuit  yields  the  differential 
relation 


V(x)  - V(x+dx)  = -dV ( x)  = [I(x)Zser  + I3Zser/2  - jcoM^ydx 
■ n<x)Zser  - I3(0-M3rr  - Zser/2)]dx 
The  corresponding  differential  equation  is 
dV(x)/dx  = -I(x)Z$er  + I3(0“M3rr  - Z$er/2) 


(6.20) 


(6.21) 


Equations  6.19  and  6.21  are  the  same  as  Eqn.  6.5,  except  that  the  term 
jwMsrr  in  Eqn.  6»5  has  been  replaced  with  the  term  (j“M3rr  + Zsgr/2)  in  the 
case  where  the  signal  rail  is  adjacent  to  the  third  rail,  and  by  the  term 
(juMarr  " zs er^  in  the  case  where  the  dc  return  rail  1s  adjacent  to  the 
third  rail.  These  changes  carry  clear  through  an  analysis  similar  to  that 
covered  in  Eqn's  6.4-6.17  above,  to  yield  these  results  for  induced  voltage  in 
the  lumped-element  equivalent  circuit  for  single-rail  track  circuits: 

V1nd  = U«M3rr  t Zser/2)I3s1nh(Td)/T  (6.22) 

where  “+"  is  for  the  case  where  the  signal  rail  is  adjacent  to  the  third  rail, 
and  is  for  the  case  where  the  dc  return  rail  is  adjacent  to  the  third 
rail . 
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6.6  Overall  Result  for  Coupling  from  Third  Rail 


The  following  relation  summarizes  the  effects  of  voltages  induced  in  track 
circuits  due  to  third-rail  currents  for  the  three  types  of  circuits: 

Vind  ■ <J“M3rr  + kZser/2)I3sinh(Yd)/Y 

k = 0 for  balanced  double-rail  track  circuit 

k = +1  for  single-rail  track  circuit  with  signal  rail  (6.23) 

next  to  third  rail 

k = -1  for  single-rail  track  circuit  with  dc  return  rail 
next  to  third  rail 

6.7  The  Conductive  Interference  Current  Transfer  Function 


Figure  6.7  shows  the  lumped-element  equivalent  circuit  with  arbitrary 
transmitting-end  and  receiving-end  impedances  lj  and  ZR  respectively.  Direct 
analysis  of  this  circuit  shows  that  the  current  into  the  receiver  IR  is 

h ’ I(d)  ■ VindY' 

(6.24) 

where  Y‘  = [Zj/IZj  + ZR)]/[(ZTI  |Zj)  + Z2  + (ZR||ZX)] 

The  current  transfer  function  relating  receiver  current  to  third-rail  current 
is  then 

HCOnd  = 1R/13  = I<d)/I3  = Y'(jwM3rr  + kZser/2)sinh( Yd)/Y  (6.25) 

For  later  use,  at  this  time  we  will  also  define  the  corresponding  transfer 
function  relating  induced  transmitter  current  to  I^: 

Hcond  - V!3  ■ I(0)/I3 

(6.26) 

= ([Z1/(Z1+ZT)]/[(ZTIIZ1)+Z2+(ZRIIZ1)]}(jwM3rr  + kZser/2)sinh(ird)/Y 
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A heuristic  understanding  of  the  nature  of  these  transfer  functions  can  be 
obtained  by  examining  the  behavior  of  Hcond  for  the  extreme  cases  of  very 
short  and  very  long  track  circuits.  For  the  case  of  the  single  balanced 
double-rail  track  circuit, 

Hcond  ' >M3rrd/CZT  + Zserd  + V d * 0 

(6.27) 

Hcond  * <j“M3rr/zser>£V(zo+ZR>]  « d * ’ 

Usually,  |Zj|  is  either  zero  or  on  the  order  of  an  ohm,  |ZR|  is  on  the  order 
of  an  ohm,  and  |ZQ|  is  on  the  order  of  a few  ohms.  Therefore  we  see  that  for 
this  case  a useful  pessimistic  approximation  is 

lHcond'  U“M3rT/zserl  " °’17  ' 176  <6-28> 

Hcond  has  been  measured  in  least  one  instance  by  Krempasky,  Clark  & 
Frasco  [Ref.  9],  for  a loop  of  rail  terminated  at  one  end  by  a running  rail  to 
running  rail  short  adjacent  to  insulated  joints,  and  at  the  other  end  by  an 
impedance  bond  adjacent  to  insulated  joints.  The  value  of  IH  ^1  found  for  a 
128*  length  of  track  at  2960  Hz  was  0.107,  as  opposed  to  the  theoretical  value 
calculated  here  of  0.17.  One  reason  surmised  for  the  discrepancy  is  the 
unknown  effect  of  possible  steel  reinforcing  rods  in  the  concrete  roadbed 
beneath  the  rails. 

For  the  case  of  balanced  2-rail  track  circuits  with  impedance  bonds,  Eqn's 
6.25-28  are  only  valid  for  the  case  of  a track  circuit  with  double  insulated 
joints  at  each  end,  and  no  third-rail  crossovers  in  the  middle.  As  was 
discussed  in  Section  2,  varying  arrangements  of  third  rail  and  impedance  bonds 
can  give  values  of  I Hcond I ^rom  zer0  t0  tw1ce  those  values.  For  instance,  in 
the  case  of  track  circuit  receiver  leads  attached  between  continuously  welded 
rails  at  the  location  of  a third-rail  crossover,  each  loop  will  contribute 
current  with  magnitude  0 . 17 1 3 to  the  conductive  interference  current  through 
the  receiver,  giving  I Hcond I = 0.34. 

Assuming  Z$er  * 6j<*>Mgrr,  the  corresponding  pessimistic  approximation  for 
single-rail  track  circuits  is 

lHcond^  ~ ^ 
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(6.29) 


Unlike  in  the  case  of  balanced  2-rail  track  circuits,  here  Eqn's  6.25-28 
generally  will  be  valid  directly  for  any  signal  block  in  which  there  are  no 
third  rail  crossovers. 

6.8  Sample  Calculations 


of 


Table  6.1  shows  the  results  of  sample  calculations  of  H 
track  circuit  length  for  the  following  cases: 


cond 


as  a function 


• Balanced  2-rail  audio-frequency  track  circuit  with  ZR  = lj  = 1 + jO  ft 
at  f = 3 kHz,  with  20  fl-kft  ballast  and  NYCTA  third-rail  spacing. 

• Single-rail  track  circuit  with  signal  rail  adjacent  to  third  rail  with 
ZR  = Zy  = 1 + jO  a at  f = 60  Hz,  with  20  n-kft  ballast  and  NYCTA 
third-rail  spacing. 

These  results  were  obtained  by  using  the  FORTRAN  77  program  shown  in  Appendix 
A.  The  purpose  of  showing  these  results  is  to  provide  an  order-of-magnitude 
glimpse  at  results,  and  to  provide  prospective  users  with  a practice  example. 


6.9  Impedance  of  the  Third  rail -DC  Return  Circuit  - Balanced  2-Rail  Case 


To  calculate  how  much  ac  current  will  flow  in  the  third  rail,  it  is 
necessary  to  know  the  total  impedance  of  all  track  circuits  hooked  end-to-end 
from  a train  to  the  substation,  as  pictured  in  Figure  6.8.  Figure  6.9  shows 
driving-point  impedances  comprised  of  track  circuits  plus  third  rail.  Three 
cases  are  shown:  balanced  two-rail  track  circuit,  single-rail  track  circuit 
with  signal  rail  adjacent  to  third  rail,  and  single-rail  track  circuit  with  dc 
return  rail  adjacent  to  third  rail.  An  extension  of  the  methods  used  in 
Sections  6. 1-6. 7 can  be  used  to  calculate  the  insertion  impedance  of  a single 
track  circuit  in  a manner  that  takes  into  account  finite  ballast  resistance 
and  the  resulting  distributed  nature  of  the  circuit.  As  will  be  shown  below, 
however,  if  the  simplifying  assumption  is  made  that  Zy  = ZR  = 0,  ballast 
conduction  ceases  to  be  a factor,  and  the  insertion  impedance  Zckt  of  a single 
track  circuit  is  then  approximately  equal  to  0.9  yh  per  meter  of  length. 
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TABLE  6.1 


Transfer  Function  Hcon(J  = IR/I3  for  Two  Track  Circuits 


Case  1 

Input  data:  d = 200  m,  f * 3 kHz,  k=0,  M3rr  = 0.234  y H/m, 
zser  = (2.00  + J23.94)  rafi/m,  Ltk  = 1.27  yH/m,  Zy  = (1  + jO)  ft  , 

ZR  = (1  + jO)  ft,  Ballast  * 20  kft-ft,  Y$h  = (164.1  + jO)  ymhos/m. 

Output  data: 

Yd  = (0.2922  + j0.2688) 

ZQ  = (8.904  + J8.192)  ft 

Zx  = (61.01  + j0.7976)  n 
Z2  = (0.2752  + J4.808)  ft 
HCOnd  = (0.1467  + j0.07321)  - 0.1640/26.52° 


Case  2 

Input  data:  d = 200  m,  f * 60  hz,  k = +1,  M3rr  = 0.234  yH/m, 

Zser  ' (190  + j692)  Ltk  " 1,67  yH/m»  ZT  " (1  + n . 

ZR  * (1  + jO)  ft.  Ballast  a 20  kft-ft,  Y$h  * (164.1  + jO)  ymhos/m. 

Output  data; 

yd  = (0.05274  + jO. 03918) 

ZQ  * (1.607  + j 1.194)  ft 

Zx  * (60.95  + jO. 02098)  ft 

Z2  * (0.03792  + JO. 1260 ) ft 

Hcond  = (0.01173  + JO. 03882)  = 0.04055/73.18° 
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DC  RETURN  CIRCUIT  WITH 


FIGURE  6.8.  Series  impedance  model  of  track  circuits  and 
third  rail  in  the  substation-to-train  loop. 
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I ( x ) - 1 3 /2 


Z 


c k t 


V!3 


a)  Balanced  double-rail  track  circuit 
I(x)-I3 

x = 0 x = d 


b)  With  signal  rail  adjacent  to  third  rail 


Z 


ckt 


I(x) 

x=0  * x=d 


With  dc  return  rail  adjacent  to  third  rail. 


FIGURE  6.9.  Series  impedance  of  one  track  circuit  and 
adjacent  thi rd  rail. 
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To  find  Z 


ckt 


= Vg/I^  as  shown  in  Figure  6.9a,  we  solve  for  I ( x)  as  a 


function  of  I3,  and  then  find  V3  by  summing  IZ  drops  in  the  third  rail,  the 
nearest  running  rail,  lower  halves  of  Z-j-  and  ZR,  and  the  induced  voltage 
caused  by  magnetic  flux  passing  between  second  and  third  rails  due  to  currents 
in  the  rails: 


V3  = CKO)  + I3/2]  ( Z-J-/2 ) + [1(d)  + I3/2]  ( Zr/2  ) 


+ V^int  + (j“V'2’')ln(d23/r3)] 


+ [Zint  + (jwy0/2Tr)ln(d23/r3)]j  [ I3/2  + I(x)]dx 

d 0 


(6.30) 


+ [(jwn0/2ir)ln(d13/d12)]' 


[ I3/2  - I(x)]dx 


In  the  above  expression,  Z3l.nt  is  the  internal  impedance  per  meter  of  the 
third  rail  and  r3  is  the  effective  radius  of  the  third  rail.  Recall  from 
Section  4 that  Zint  was  the  internal  impedance  per  meter  of  the  running  rail 
and  r^  and  r2  were  the  effective  radii  of  the  running  rails.  If  the 
techniques  described  in  Section  4 have  been  used  to  measure  rail  impedance, 
then  r1  9 - are  taken  to  be  equal  to  the  radius  of  the  reference  pipe,  and  the 
internal  impedances  are  equal  to  the  measured  impedances. 


Integration  of  both  sides  of  Eqn.  6.5  yields  the  relation 
d 

I(x)dx  = (l/Z$er)[V(0)  - V(d)  + j“M3rrdI3] 

) 

Then,  since 


v(o)  = -i(o)zT  = -iTzT 

V(d)  = I(d)ZR  = IRZR  = 


■ _I3HcondZT 
!3HcondZR 


and 


where  Hcond  and  H^ond  are  taken  from  Eqns  6.25,26  with  k = 0, 
the  insertion  impedance  Zckt  of  the  track  circuit  reduces  to: 


(6.31) 


(6.32) 


-50- 


Zckt  ■ ‘W 


■ (ZT  + ZR>/4  + dlZ3int  + <j“%/2’)ln(d23/r3) 


+ Zint/2  + (jwM0/4ir)ln(d23/r2)  + ( jaiyo/4Tr)l  n(d13/d12) } 


+ H'  1J2  + H ,Zn/2 
cond  T cond  R 


(6.33) 


((J“M3rrd  * 


U - 


[Zint  + (j“M0/2ir)ln(d23/r2)  - (jwy0/2ir)ln(d13/d12)]} 

To  gain  a bit  of  insight  into  the  magnitude  of  Zckt,  it  is  useful  to  look 
at  a simplifying  case.  If  we  set  Zj  = ZR  = 0,  and  assume  that 

r3  = r2*  Z3int  = Zint*  and  ^“M3r/Zser^  = 1^6,  the  above  expression  can  be 
written  in  terms  of  an  impedance  per  meter: 


Wd  • “! 5/3)Zint 

(6.34) 

+ (j»u0/2.)[(5/3)ln(d23/r2)  + ( 1/3 ) 1 n( d13/d12 ) D } 

If  we  then  define  Zser23  to  be  the  series  transmission  line  impedance  per 
meter  of  the  transmission  line  composed  of  the  second  and  third  rails,  the 
Zint  and  1°9  terms  in  the  above  expression  are  replaced  by  a Zser03 

term,  yielding  the  expression 


Wd  ■ (5/6)Zser23  + 


( 1/3  )fj  <oyo/2i)l  n ( d 


13/d12) 


(6.35) 


Assuming  for  the  moment  that  Z^  = 0,  and  that  rail  diameter  and  spacing 
are  typical  NYCTA  values,  namely  d23  = 67.3  cm,  d^2  = 150  cm,  d^3  = 217  cm, 
and  r2  = 5 cm,  the  resulting  value  of  l_ckt  = Zckt/jwd  is  0.9  microhenries  per 
meter. 
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6.10.  Impedance  of  the  Third  Rail -DC  Return  Circuit  - Signal  Rail 
Adjacent  to  Third  Rail 


Figure  6.9b  shows  a pair  of  running  rails  with  a single-rail  track  circuit 
in  which  the  signal  rail  is  adjacent  to  the  third  rail.  In  order  to  obtain 
the  voltage  V3  due  to  current  I3,  we  sum  voltages  along  path  ABCDEF  and  obtain 
the  relation 


V3  = I(0)ZT  + 1(d) ZR 


+ I3d*{Z3int.  + (ju>M0/2ir)[ln(d23/r3)  + ln(d13/d12)]} 


(6.36) 


+ {Z  + (j«u  /2ir)[ln(d  Jr  ) -ln(d  /d  )]} 
int  0 23  2 13  12 


I(x)dx 


We  now  integrate  both  sides  of  Eqn.  6.19  to  obtain  the  relation 
d 

fl(x)dx  = J[V(0)  - V(d)]/Zserl  + I3d-(j^3rr/Zser  + 1/2) 

0 

Then  using  Eqn's  6.25,26  with  k = +1/2,  we  write 


(6.37) 


V(0)  - -I(0)ZT  - -H'ondZTI3 


V(d)  - +I(d)ZR  - HcondZRI3 


(6.38) 


Substituting  the  results  of  Eqn.  6.38  into  Eqn.  6.37,  and  the  result  of 
that  into  Eqn.  6.36,  we  obtain  the  following  expression  for  the  insertion 
impedance  of  the  track  circuit  for  this  case: 


zckt  ' <W 


(6.39) 


- d*lZ3int  + (j“M0/2»)[ln(d23/r3)  + ln(d13/d12)]) 


+ d'(j“M3rr/Zser  + + ( j »>u0/2ir ) [ 1 n( d23/r2 ) - ln(d23/d12]) 

+ (HcondZT  + HcondZRU1  ' (1/Zser)[Zint  + « j“V2’>ln(d23/r2> 


- (j^y0/2ir)ln(d13/d12)]} 


-52- 


It  so  happens  that  for  the  case  where  Zj  = ZR  = 0,  Eqn's  6.33  and  6.39  can 
be  put  in  identical  form,  and  the  track  circuit's  impedance  characteristics 
can  once  again  be  described  by  an  approximate  inductance  per  unit  length: 
Lckt  = Zckt^a)d  ~ 0,9  m1’crohenries  per  meter. 

6.11.  Impedance  of  the  Third  Rail -DC  Return  Circuit  - DC  Return  Rail 
Adjacent  to  Third  Rail 


Figure  6.9c  shows  a pair  of  running  rails  with  a single-rail  track  circuit 
in  which  the  dc  return  rail  is  adjacent  to  the  third  rail.  In  order  to  obtain 
the  voltage  V3  due  to  current  I3  in  this  case,  we  sum  voltages  along  path  ABC 
and  obtain  the  relation 


v3  = V^int  + Z1nt>  + U“V2l')[ln(d23/r3)  + ln(d23/r2)]) 


d 

» 


(6.40) 


{Z.  . + 
1 int 


(juu0/2ir)[ln(d23/r2)  - ln(d13/d12)]}* 


I ( x ) dx 


0 

We  now  integrate  both  sides  of  Eqn.  6.21  and  for  this  case  obtain  the 
relation 
d 


I ( x) dx  = tCVCO)  - V(d)]/Zser  + I3d-(j“M3rr  - Zser/2)/Zser  (6‘41) 

4 

0 


Then  using  Eqn's  6.25,26  with  k = -1/2,  we  write 


V(0)  = -I(0)ZT  = -H‘ondZTI3 
V ( d ) - +1 C d)ZR  = HcondZRI3 


(6.42) 


Substituting  the  results  of  Eqn.  6.42  into  Eqn.  6.41,  and  the  result  of 
that  into  Eqn.  6.40,  we  obtain  the  following  expression  for  the  insertion 
impedance  of  the  track  circuit  for  this  case: 
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= d’{Z3int  + Z1nt  + U“*‘</2,r)[1n(<l23/r3>  + Inl^/r^]) 


(6.43) 


hrr^ser  " 


{Zint  + ^'a,V2ir)Cln(d23/r2)  " ln(di3/di2)]} 


Once  again,  if  we  allow  ZR  = Zj  = 0,  Eqn.  6.43  gives  exactly  the  same 
result  as  do  Eqn's  6.33  and  6.39  when  ZR  = lj  = 0,  and  the  same  value  of  Lckt 
of  0.9  microhenries  per  meter. 

6.12  Rail  Data  to  Use  in  Computation 

The  equations  above  in  Section  6 call  for  rail  radii  and  internal 
impedances.  There  are  two  possible  ways  to  choose  the  radii  and  impedances  to 
use. 

The  first  way  to  proceed  is  to  use  the  measurements  of  rail  impedance 
outlined  in  Section  3 as  a direct  source  of  rail  impedance  data.  In  this 
case,  the  rail  radii  r^,  r^,  and  r^,  are  all  set  equal  to  rcu,  the  radius  of 
the  reference  pipe  used  in  the  measurements.  Then,  internal  rail  impedance 
Z.  ^ is  set  equal  to 


The  corresponding  expression  is  used  to  find  1 . for  the  third  rail.  Z 

ji nt  S6r 

can  be  set  directly  equal  to  Rtk  + from  Section  3. 

A second  way  to  proceed  is  to  use  the  graphical  procedure  outlined  in 
Section  5 to  find  effective  values  for  rail  radii  and  (m/o)  ratios.  These 
values  are  used  for  r^,  r2,  and  r^.  Then  Eqn.  5.6  is  used  to  calculate 
approximate  effective  values  of  internal  rail  resistance  and  inductance.  This 
second  way  of  proceeding  is  also  based  on  the  original  rail  measurement  data. 


Z.  . = Z /w  ^ R + jojAL/w 
int  meas  rail 


meas 


(6.44) 
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but  the  data  is  fitted  to  a theoretical  model  of  rail  impedance  behavior  that 
is  then  used.  This  second  method  can  be  thought  of  as  a shorthand  method, 
since  the  original  rail  data  are  condensed  to  a value  of  radius  and  (u/°). 

6.13  Summary  & Conclusions 


This  section  has  introduced  the  analytical  methods  for  calculating  the 
conductive-interference  coupling  between  third  rail  and  track  receivers,  in  a 
circuit  model  that  includes  distributed  ballast  conduction  effects.  Rather 
idealized  situations  have  been  dealt  with.  Such  factors  as  arbitrary  third- 
rail  geometry  have  been  neglected.  The  existence  of  third-rail  crossovers  at 
arbitrary  locations  can  be  dealt  with  by  joining  equivalent  n-circuits  as 
shown  in  Figure  6.6,  without  including  or  ZR. 

One  major  result  of  the  analysis  presented  in  this  section  is  the 
prediction  of  the  impedance  of  the  circuit  formed  by  running  rails  and  dc 
return  rail.  This  impedance  is  a series  element  in  the  circuit  extending  from 
substation  to  train  in  which  audio-frequency  substation  harmonic  currents  and 
chopper  harmonic  currents  flow.  The  simplified  model  that  assumes  that  lj  and 
ZR  are  both  0 for  every  track  circuit  is  also  a good  worst-case  model.  This 
model  predicts  that  the  series  inductance  of  this  circuit  is  approximately  0.9 
microhenries  per  meter  of  circuit  length--a  value  smaller  than  but  comparable 
to  the  running-rail  transmission-line  inductance  which  has  a nominal  value  of 
approximately  1.4  microhenries  per  meter. 

The  0.9  uH/meter  figure  is  a conservative  one,  and  can  therefore  be  used 
in  preliminary  analysis  to  find  conservative  estimates  of  expected  Cl. 
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7.  INTERFERENCE  SOURCE  CHARACTERISTICS  OF  MULTI-CAR  TRAINS:  CIRCUIT  EFFECTS 


7.1  The  Multi -Car  Source 

The  following  sections  describe  the  characteristics  of  multi -car  rapid 
transit  trains  as  sources  of  conductive  interference.  Experience  gained  to 
date  has  been  with  cars  powered  by  chopper-controlled  dc  propulsion  systems. 
However,  the  analytical  techniques  described  are  directly  applicable  to 
solid-state  controlled  ac  propulsion  systems.  This  section  discusses  the 
circuit  aspects  of  multi -car  trains.  The  following  section  discusses 
statistical  characteristics. 

A multi -car  train  consists  of  a number  of  practically  identical  sources  of 
conductive  interference  generating  signals  that  add  together.  Each  car  in  a 
train  will  have  very  nearly  the  same  chopper  frequency  as  the  other  cars, 
differing  only  by  tolerances  of  the  frequency  standards.  The  chopper  harmonic 
components  of  each  car  will  have  approximately  the  same  amplitude  vs.  time 
characteristics  as  the  other  cars  in  the  train. 

A multi -car  train  is  a distributed  source  having  the  form  of  a ladder 
network.  Each  car  provides  a rung  consisting  of  the  car's  Norton  equivalent 
chopper  harmonic  output  current  source  in  parallel  with  a source  inductance. 
Inter-car  inductance  provides  the  series  element  in  the  circuit  between  each 
rung.  Since  the  interference  current  from  any  one  car  is  attenuated  as  it 
travels  down  the  ladder,  cars  will  contribute  slightly  different  magnitudes  of 
current  to  the  third  rail  segment  running  from  the  end  of  the  train  to  the 
substation.  Truman  [Ref.  10]  has  used  this  model  to  calculate  the  circuit 
characteristics  of  multi -car  trains. 

7.2  The  Equivalent  Circuit  of  a Multi -Car  Train 


Figure  7.1  shows  a model  for  the  distributed  circuit  of  a multi -car  train 
and  a third-rail  running-rail  loop  ending  at  a substation.  Each  car  is 
represented  by  its  Norton-equivalent  current  source  and  source  inductance. 
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FIGURE  7.1.  The  equivalent  circuit  for  conductive 

interference  generation  in  a multi-car 
train. 
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The  source  inductance  is  assumed  to  be  equal  to  the  car's  line  filter 
inductance.  The  series  inductance  between  each  train  and  its  neighbors  is 
assumed  to  be  equal  to  the  series  inductance  per  meter  of  the  third-rail 
running-rail  loop  times  the  car  length. 

This  model  uses  simplifying  assumptions  that  are  not  believed  to  harm  the 
qualitative  or  quantitative  nature  of  the  results.  In  reality,  each  car  is 
attached  to  the  third  rail  by  two  shoes,  and  to  the  running  rails  by  8 wheels. 
Furthermore,  the  series  inductance  between  cars  in  the  circuit  will  depend  on 
the  placement  of  propulsion  conductors  within  the  car. 

In  Figure  7.1,  the  car  produces  current  I .,  and  has  car  output 

C J 

inductance  L . The  series  inductance  between  one  car  and  the  next  is  L . 

car  cc 

The  total  Norton  equivalent  source  current  of  an  n-car  train  is  In,  and  the 

Norton  equivalent  source  inductance  of  the  n-car  train  is  Ln.  The  values  of 

L for  n = 1,  2,  . . . can  be  found  by  the  recursion  relation 
n 

^n  “ ^n-1  + *"cc^  ^car*  *"1  “ ^car  (7.1) 

A similar  recursion  relation  can  be  used  to  find  the  expression  for  source 
current: 


I = I + I . L ./(l  , + L ), 
n cn  n-1  n-1  n-1  cc 


= ‘cl 


^n  anl^cl  + an2*c2  + * * * + ann^cn*  ann  * 


(7.2) 


By  way  of  example,  these  relations  were  used  to  calculate  Ln  and  In  for 
values  of  n between  1 and  8,  assuming  L.,„  = 500  yh  and  L = 20  uh: 

Luff  LL 
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L1=500  yh,  I1=Icl 

L2=255  yh,  I2=,962Icl+Ic2 

L3=177.4  Mh,  I3=.892Ic1+.928Ic2+Ic3 
L4=141.5  nh,  I4=.801Ic1+.882Ic2+.898Ic3+Ic4 

(7.3) 

L5-122.1  Ph,  I5=.701Ic1+.773Ic2+.787Ic3+.876Ic4+Ic5 
L6=110.6  Ph,  I6=.602Ic1+.664Ic2+.676Ic3+.753Ic4+.859Ic5+Ic6 
L?=103.6  Ph,  I7=.510Ic1+.562Ic2+.572Ic3+.638Ic4+.727Ic5+.847Ic6+Ic7 

l8=  "-1  “h-  I8“-427Icl+-471Ic2+-497Ic3+-535Ic4+-6091c5+-710Ic6+-838Ic7+Ic8 

In  the  above  set  of  relations,  Ln  appears  to  be  approaching  an  asmyptotic 
value  for  large  n.  This  value,  defined  as  L^,  can  be  found  by  solving  the 
relation 


L®  ^L°°  +Lcc)  * *Lcar 


(7.4) 


For  L = 20  yh  and  = 500  yh,  the  resulting  value  is  L = 90.5  yh. 
cc  car  °° 

We  will  assume  that  all  the  Icn's  are  phasors  of  equal  amplitude  Ic>  and 
random  phase  with  respect  to  each  other.  The  maximum  magnitude  of  In,  defined 
as  Inmax  will  occur  when  all  of  the  individual  contributions  add  in  phase: 


^ax  ~ Vd 


a ) 
cn 


For  the  coefficients  given  in  Eqn.  7.3, 


(7.5) 


lmax  " *c 

L av  = 1.9621 
2max  c 

^ax  ' 2‘82Ic 

!4max  ■ 3*581 

5m ax  = 4*14Ic 

r6max  = 4,55Ic 

*7max  ‘ 4'86Ic 

*8max  ' 5'071 
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As  n grows  large,  the  value  of  I also  approaches  an  asmyptotic  value 

I II II 01 A 

v»  found  by  solving  the  relation 

“II  a a 

^ax  *c  + + Lcc)  (7.7) 

to  obtain  the  result 


I 


«max 


*c^«  + Lcc^Lcc 


(7.8) 


For  the  values  = 90.5  yh  and  Lcc  = 20  yh,  the  resulting  value  is 
I = 5.52I„.  This  asmyptotic  behavior  results  when  each  car  added 

“fllaX  C 

attenuates  as  much  of  the  maximum  current  of  the  other  cars  as  it  then 
contributes  back. 


It  is  also  instructive  to  calculate  the  rms  value  of  I , noted  as  I 

n nrms 

Since  the  individual  components  of  In  have  random  phase  with  respect  to  each 
other,  the  mean  squared  and  rms  values  of  In  are 


I2  = I2(a2  + . . . + a2  ) 

nrms  cv  nl  nn^ 


I _ t /p2  , , 2 x 1/2 

nrms  c nl  nn 


(7.9) 


Based  on  the  coefficients  from  Eqn.  7.3  in  the  example  above. 


*lrms 

‘c 

^nms 

■ 1-37IC 

*3rms 

■ !-63Ic 

’4ms  ■ U801c 

(7.10) 

I_ 

1.861 

I _ 

= 1.891 

I_ 

= 1.891 

ln  = 1.871 

5rms 

c 

6rms 

c 

7 rms 

c 

8 rms  c 

In 

the  data 

above. 

I__e  is  seen  to 
nrms 

peak  and 

then  start  declining 

toward 

the  asmyptote  1^  . An  expression  for  I is  found  by  squaring  each  term 

in  an  expression  equivalent  to  Eqn.  7.7  to  obtain  the  relation 


I 


2 

00  rms 


rc  + I«rms[L-/(L»  + Lcc)]; 


(7.11) 
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This  relation  has  the  solution 


= ’r[l  - l£/CL_  + L )2]1/2  (7.12) 

®rms  c oo  oo  cc 

For  the  numerical  data  used  in  the  example  above,  the  surprisingly  small 

result  is  r e = 1.741  . 

°°rms  c 


To  find  the  actual  third-rail  current  Ig  that  flows  at  a given  chopper 
harmonic  frequency,  one  must  consider  the  circuit  shown  in  Figure  7.2  that 
includes  the  Norton  equivalent  source  for  the  train,  and  the  third-rail 
running-rail  circuit  to  the  substation.  The  Norton  equivalent  source  current 
In  divides  between  the  train's  source  inductance  Lp  and  circuit  inductance 
L3rrD*  The  resu1tin9  third-rail  current  is 


*3  = W(Ln  + L3rrD> 


■ VC1  + <L3rrD'Ln>} 


(7.13) 


Examination  of  this  equation  shows  that  as  the  number  of  cars  in  the  train 
grow  larger  and  Ln  consequently  grows  smaller,  the  ratio  Ig/In  grows  smaller. 
This  behavior  is  due  directly  to  the  fact  that  each  car  in  the  train  provides 
an  alternate  shunt  path  for  the  interference  current  from  each  other  car. 
Increasing  the  value  of  LgrrD,  the  inductance  of  third-rail  running-rail 
circuit,  also  decreases  the  ratio  I3/ In . 

One  example  is  given  here  of  a value  of  third-rail  current  Ig  for  a multi - 
car  train.  We  use  the  sample  data  given  above  with  Lg  = 99.1  yh;  we  use  the 
value  Lgrr  = 0.9  yh/meter  found  in  Sec.  6;  and  we  choose  D = 100  meters.  Then 
the  circuit  constant  from  Eqn.  7.13  becomes  Lg/(Lg  + LgrrD)  = 0.524;  and  solv- 

in9  for  '3raax  by  lettin9  h = W in  E9n-  7-13-  we  obtain  !3nlax  = 2-66IC 

Correspondingly  smaller  values  of  Inmax  are  found  for  smaller  values  of  n, 
and  for  larger  values  of  D. 
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FIGURE  7.2.  The  conductive  interference  circuit 
showing  train,  circuit  inductance, 
and  substation. 
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7.3  Summary  and  Conclusions 


In  conclusion,  we  note  that  a multi -car  train  has  a significantly  smaller 
source  inductance  than  a one-car  train.  Therefore,  the  addition  of  a very 
moderate  amount  of  circuit  inductance--even  as  little  as  that  of  100  meters  of 
track— can  significantly  reduce  maximum  values  of  chopper  harmonic  current. 
This  fact  should  be  investigated  as  a possible  means  of  limiting  conductive 
interference  current  at  those  locations  where  intolerably  large  values  are 
encountered  when  trains  operate  near  substation  feeds. 

Power  circuit  designers  generally  try  to  minimize  the  inductance  of 
circuits  carrying  large  magnitudes  of  dc  current,  because  high  circuit 
inductance  increases  the  difficulty  of  interrupting  fault  currents  in  dc 
circuits.  However,  the  amount  of  added  circuit  inductance  recommended  here  is 
very  small  compared  to  the  overall  inductance  in  long  third-rail  circuits. 

The  analysis  in  this  section  only  deals  with  levels  of  conductive 
interference  contributed  by  choppers  on  the  rapid  transit  cars.  Substation 
conductive  interference  is  more  directly  dealt  with.  Given  the  levels  of 
circuit  inductance  in  question,  rectifier  substations  are  nearly  ideal  voltage 
sources.  As  n increases  and  Ln  decreases,  the  total  inductance  of  the  circuit 
in  series  with  the  substation  decreases,  and  the  level  of  I3  at  substation 
harmonics  increases  accordingly. 
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8.  INTERFERENCE  SOURCE  CHARACTERISTICS  OF  MULTI-CAR  TRAINS:  STATISTICS 


8.1  Statistical  properties 


Each  car's  chopper  control  circuitry  is  independent  of  other  cars. 
Therefore,  at  each  chopper  harmonic  frequency  in  question,  the  conductive 
interference  signals  from  the  different  cars  of  a train  will  be  randomly 
phased  with  respect  to  each  other.  In  a long-term  statistical  sense,  at  any 
chopper  harmonic  frequency,  the  total  average  noise  power  will  equal  the  sum 
of  noise  powers  from  the  individual  cars.  In  other  words,  the  total  rms 
current  will  equal  the  square  root  of  the  sum  of  the  squares  of  the  individual 
components.  However,  at  some  time  the  individual  components  can  be  expected 
to  add  practically  in  phase  to  produce  a signal  whose  amplitude  is  the  direct 
sum  of  the  individual  components.  There  will  be  a statistical  distribution  of 
interference  levels  characterized  by  a probability  density  function  taking  on 
non-zero  values  from  zero  to  the  maximum  interference  level. 


As  a function  of  time,  a train's  Norton-equivalent  conductive  interference 
source  current  at  a particular  chopper  harmonic  frequency  f = <*>o/2ir  can  be 
written  as 


'n141  = Ic[anlcos(“ot  + 61>  + 


+ acos(w  t + el] 
nn  o n 


= I r cos(u  t + e) 
c n o 


(8.1) 


As  noted  here,  the  total  conductive  interference  source  current  for  an 
n-car  train  is  equal  to  the  sum  of  n independent  phasor  quantities.  The 
frequency  of  any  particular  chopper  harmonic  typically  will  vary  by  a fraction 
of  a hertz  from  car  to  car.  This  will  cause  the  values  of  e^,  e^ , ...»  in  the 
equation  above  to  vary  slowly  with  time.  Therefore,  over  the  short  time  span 
it  takes  for  a train  to  accelerate  from  rest  up  to  the  speed  at  which  the 
choppers  stop  chopping,  phase  relations  between  the  interference  current 
phasors  of  cars  at  a particular  chopper  harmonic  will  remain  approximately 
fixed.  However,  over  the  longer  period  between  starts,  phases  will  have  a 
chance  to  re-randomize. 
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The  statistics  of  the  magnitude  of  I depend  on  the  statistics  of  the 
phasor  additions  of  the  terms  anj.cos(uQt  + 9j)  to  produce  the  phasor  whose 
magnitude  is  r . The  computational  problem  has  been  simplified  by  considering 
the  single  limiting  case  in  which  all  the  a ,'s  have  unity  value.  This  case 

* I J 

is  equivalent  to  the  case  in  which  the  intercar  inductance  L is  zero. 

V V 

8.2  Probabil ity  , Density  Functions  [Ref.  11] 

Figure  8.1a  depicts  the  addition  of  n randomly-oriented  unit-amplitude 

phasors  for  the  case  n = 4.  Figure  8.1b  shows  the  phasor  addition  of  a phasor 

of  unit  amplitude  to  a phasor  of  arbitrary  amplitude  r'  = rp  ^ to  produce  a 

new  phasor  of  amplitude  r = rn , when  the  phase  difference  between  r'  and  the 

unit  phasor  is  e.  For  conciseness,  we  use  r in  place  of  r , and  r'  in  place 

of  r , in  the  equations  below. 
n-1 

To  develop  the  probability  density  functions  for  r , an  iterative 
technique  was  developed  that  allows  the  probability  density  function  pn(rn)  to 
be  calculated  from  the  probability  density  function  pn_1 ( rn_1 ) = pn_1 ( r ' ) . 
Note  in  Fig.  8.1b  that 

r = [ ( r ' + cose)2  + sin2e]1^2  (8.2) 

2 2 

Defining  c = cose  and  use  of  the  relation  cos  e + sin  e = 1 allows  us  to  write 

r = [1  + r'2  + 2r'c]1//r2  (8.3) 

This  relation  can  be  transformed  into  one  giving  c in  terms  of  r*  and  r: 

c(r,r')  = (r2  - 1 - r'2)/2r‘,  for  r,r'  > 0 (8.4) 

For  any  value  of  r > 0,  the  probability  density  function  p (r)  is  defined  by 
the  following  expression,  which  states  that  the  probability  that  a sample 
value  of  r,  denoted  by  R,  lies  between  r and  r+dr  is  equal  to  pn ( r ) dr : 

P[r  < R - r+dr]  = pn(r)dr  (8.5) 
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a)  Addition  of  4 randomly  oriented  unit  phasors. 


b)  Addition  of  a unit  phasor  of  random  orientation  to 
a phasor  of  random  length  to  produce  phasor  r. 


FIGURE  8.1.  Addition  of  randomly  oriented  phasors. 
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Regarding  c and  r'  as  statistically  independent  random  variables,  R will  have 
a value  between  r and  r+dr  provided  r'  and  c jointly  have  values  in  a strip  in 
the  (r',c)  plane,  as  shown  in  Figure  8.2.  The  lower  edge  of  the  strip  is 
given  by  c = c(r,r'),  and  the  upper  edge  by  the  relation  c+dc  = c(r+dr,r'). 
The  differential  dc  is 


dc  = (3c/ar)dr  = (r/r')dr 


(8.6) 


It  follows  that 


Pn(r)dr  = P[r  < R - r+dr]  = 


Pn_i ( r ' ) Pc ( c ) dr ' dc 
over  strip 


(8.7) 


dr*  »pn _1(r' ) pc (c(r,r* )) dc 
r'=0 


The  parameter  c has  statistics  that  obey  the  inverse-cosine  probability 
law.  This  law  gives  the  probabilities  of  results  when  the  function  c = cos(e) 
is  sampled  a large  number  of  times  for  statistically  independent  values  of  0: 


P ( C ) = l/ir(  1 - C2)1/2,  for  -1  < c < +1 

V 


= 0, 


otherwi se 


(8.8) 


When  the  form  for  c given  by  Eqn.  8.4,  the  functional  form  of  p ( c ) from  Eqn. 

V 

8.8,  and  the  form  for  dc  given  by  Eqn.  8.6  are  substituted  into  Eqn.  8.7,  the 
resulting  relation  for  the  pn(r)  becomes 

1+r 

Pn(r)  = (2r/ir)  pn_1  ( r ' ) [4r ' 2 - ( r2-r ' 2-l)2]”1^2dr'  (8.9) 

ll-rf 

The  limits  in  the  above  integral  are  found  by  setting  c = -1  in  Eqn.  8.4  and 
solving  for  r'  in  terms  of  r. 
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c ( r + d r , r ' ) 
c (r  ,r'  ) 


FIGURE  8.2. 


The  strip  in  the  (r',c)  plane  where 
r < R £ r+dr. 
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This  relation  can  be  applied  immediately  to  the  known  form  for  n = 1, 
which  is 

p^r)  = 6 ( r- 1 ) (8.10) 

to  yield  the  result  for  n = 2.  Substituting  p^r')  = <5(r'-l)  into  Eqn.  8.9 
and  integrating,  we  obtain 

P2(r)  = 2Ar<4  - r2)1/2, 

= 0,  otherwise 

To  continue  beyond  the  case  n = 2 requires  the  use  of  a digital  computer 
to  perform  numerical  integration.  Changes  of  variables  were  made  to  remove 
singularities  in  functions  under  the  integral  sign,  and  other  analytical  and 
numerical  techniques  were  used  to  improve  accuracy.*  Then,  the  Pn(r)  func- 
tions were  found  numerically  for  n = 3,..., 8.  The  resulting  Pn(r)  functions 
for  n = 1,...,8  are  shown  in  Figure  8.3.  Hoelscher  & Rudich  have  obtained 
equivalent  plots  for  p^r)  - pQ(r)  by  Monte  Carlo  techniques  [Ref.  8]. 

8.3  Statistical  Results 


for  0 - r < 2 


(8.11) 


Theoretical  and  analytical  results  for  the  statistics  of  rp  are  summarized 
below: 

p (r  ) = 0,  for  r < 0 and  for  r > n 
n n n n 

n 

> 

pn(rn)drn  = l>  for  a11  n 
')  n 

(Mean  squared  value  of  rR)  = <r2>  = jrnPn(rn)drn  = n 

0 

(RMS  value  of  rn)  = <r^>1^2  = n1/2  (8.12) 

n 

(Average  value  of  r ) = <r  > = r p (r  )dr 

n n J n n n n 

0 

* The  author  indebted  to  Tatiana  Vinnikova,  for  developing  changes  in  vari- 
ables that  removed  singularities  in  the  integrand  of  Eqn.  (8.9)  at  the 
limits  of  integration,  and  for  performing  the  numerical  integrations. 
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r 


FIGURE  8.3.  Probability  density  functions  pk  ( r ) for  sums  of  k 

randomly  oriented  unit  phasors,for  k = 1,  ...»  8. 
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(8.12,  cont'd) 


<r^>  = 1.80 
<r7>  = 2.37 


<r2>  = 4/tt  = 1.27 
<rg>  = 2.01 
<rg>  =2.53 


8.4  Behavior  for  Large  Values  of  n 


<rg>  = 1.58 
<rg>  = 2.20 


For  n » 1,  the  probability  density  function  p ( r ) approaches  the  Rayleigh 

2 n 

probability  density  function,  with  parameter  a = n/2: 


PR(r)  = (r/a2)exp(-r2/2o2) , for  r - 0 
= 0,  for  r < 0 


(8.13) 


In  other  words,  as  n grows  large,  the  sum  of  n randomly  phased  unit-amplitude 
phasors  has  an  envelope  whose  statistics  are  the  same  as  the  envelope  function 
of  a narrow-band  gaussian  random  process.  This  fact  is  shown  graphically  in 

p 

Fig.  8.4,  in  which  pQ( r ) and  the  Rayleigh  pdf  pR(r)  = (r/4)exp(-r  /8)  are 
plotted  together. 

8.5  Probability  of  r Exceeding  a Particular  Value 


The  function  that  gives  the  probability  that  a sample  value  R is  less  than 


or  equal  to  r is  the  probability  distribution  function  Pn( r ) : 


P[R  - r]  = Pn(r)  = 


pn(x)dx 


(8.14) 


We  can  take  the  lower  limit  of  the  above  integration  to  be  zero. 


The  function  that  gives  the  probability  that  a sample  value  R is  greater 


than  r is  the  complementary  probability  distribution  function  P^(r): 


P[R  > r]  = P^(r)  = 


pn(x)dx 

% 

Numerically  derived  plots  of  P^( r ) are  shown  in  Fig.  8.5  for  n = 2,..., 8. 


(8.15) 
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RAYLEIGH 


r 


FIGURE  8.4.  The  Rayleigh  pdf  with  a2  = 4 compared 
to  p8(r). 
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r 


FIGURE  8.5.  Complementary  probability  distribution 

functions  P ^ ( r ) for  sums  of  k unit  phasors. 
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8.6  Extension  to  the  Case  of  Unequal  Contributions  From  n Cars 


As  the  results  of  Section  7 showed,  the  cars  of  an  n-car  train  all 
contribute  different  amounts  to  the  total  conductive  interference.  We  may 
regard  the  weighting  functions  for  the  various  cars  in  an  n-car  train  to  be 
phasors  with  amplitudes  anl,  ...  , ann,  and  statistically  independent  phases. 
Then,  if  we  really  want  to  find  pn(r)  for  a particular  value  of  n,  here  is 
what  we  do: 

Starting  with  the  first  phasor  with  amplitude  aRl  and  pdf 
px(r)  = 6(r  - anl) 

we  recursively  use  the  following  integral  relation  to  find  succeeding  pdf's  up 


Equation  8.16  was  obtained  by  a procedure  directly  analagous  to  that  followed 

in  Sec.  8.2.  A given  probability  function  Pn(r)  will  have  zero  value  for  r<0 

and  for  r > (a  . + ...  + ann) . 

nl  nn 

Note  that  unlike  the  case  where  all  the  a.  .'s  have  unit  magnitude,  we 

0 

cannot  use  the  pdf  for  an  n-car  train,  to  find  the  final  pdf  for  an  (n+l)-car 

train.  We  have  to  start  from  scratch  with  the  set  of  a.  coefficients  for 

jn 

each  value  of  n,  and  work  our  way  through  intermediate  pdf's  to  find  the  final 

form  of  p ( r) . 
n 

2 2 2 

We  do  know  that  for  any  particular  value  of  n,  <rn>  = (anl  + ...  + ann). 

We  can  surmise  that  as  n grows  large,  pn(r)  approaches  a Rayleigh  pdf  with 

2 2 

parameter  a = <r  >/2. 


to  pn(r): 


r'2  - a2  )2]-1^2dr 
nk 


(8.16) 
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8.7  Statistics  Summary 


The  statistical  results  above  represent  formal  answers  to  a highly 
idealized  problem.  Nonetheless,  these  results  can  serve  as  a guide  to  our 
thinking  about  expected  conductive  chopper  interference  levels  from  multi -car 
trains.  Clearly,  the  writers  of  specifications  for  limits  on  conductive 
chopper  interference  must  be  precise  in  the  statement  of  what  their  limits 
pertain  to.  For  instance,  when  limits  for  one  car  are  extrapolated  to  an 
n-car  train,  we  must  state  whether  the  the  limit  for  the  train  is  a maximum 
figure  or  an  rms  figure.  In  addition,  we  must  remember  that  in  this  problem, 
"nms"  refers  to  an  averaging  process  over  a statistical  ensemble,  and  not  to 
an  averaging  process  over  the  period  of  a predictable  periodic  signal. 

We  saw  in  Section  7 that  at  a particular  chopper  harmonic,  the  maximum 
conductive  interference  current  flowing  in  the  third  rail  from  an  n-car  train 
would  be  a fraction  of  the  current  from  one  car  times  n.  In  this  section,  we 
saw  that  actual  values  of  third-rail  current  will  be  randomly  distributed 
between  zero  and  that  maximum  value. 

Perhaps  it  is  of  comfort  when  measuring  conductive  chopper  interference 
from  multi -car  trains  in  the  field,  to  realize  that  the  randomness  in  the  data 
is  not  necessarily  due  to  experimental  error,  but  is  at  the  heart  of  the 
process  itself. 
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9.  MEASUREMENT  OF  CONDUCTIVE  INTERFERENCE:  TECHNIQUES  & DATA 


9.1  Recommended  Practices 


A draft  set  of  recommended  practices  has  been  developed  under  the  aegis  of 
the  Rail  Transit  Electromagnetic  Interference  Technical  Working  Group  (TWG). 
This  committee's  membership  includes  representatives  of  rapid  transit 
signaling  and  propulsion  manufacturers,  transit  industry  consultants,  and 
university  and  government  researchers.  TWG  is  organized  under  the  auspices  of 
DOT/UMTA.  TWG  has  prepared  draft  Recommended  Practices  that  describe  methods 
for  observing,  recording,  and  analyzing  conductive  interference  in  rapid 
transit  signaling  systems,  and  for  measuring  the  susceptibility  of  signaling 
systems  to  conductive  interference  [Ref.  12].  The  techniques  described  below 
for  observing  and  recording  conductive  interference  are  described  in  large 
part  in  the  Recommended  Practices.  Measurement  techniques  are  summarized 
below. 

9.2  Instrumentation 


Conductive  interference  in  rapid  transit  signaling  is  in  the  form  of 
power-frequency  or  audio-frequency  current  components  superimposed  on  dc 
propulsion  current.  Propulsion  current  levels  can  be  large--thousands  of 
amperes.  Signaling  current  levels  are  small --amperes  in  the  case  of 
power-frequency  signaling,  or  tens  to  hundreds  of  milliamperes  in  the  case  of 
audio-frequency  signaling.  Therefore,  special  means  must  be  used  to  observe 
conductive  interference  currents  whose  magnitudes  will  be  of  the  same  order  as 
those  of  signaling  currents. 

The  primary  tool  used  for  measuring  conductive  interference  currents  is 
the  Hall -effect  servo  current  sensor.  The  type  used  in  this  activity  is  the 
LEM  Transfoshunt} . Transfoshunt  models  having  maximum  dc  current  capability 
of  1200  amperes  and  10,000  amperes  have  been  used  in  this  program--the  lower 
current  model  aboard  single  rapid  transit  cars,  and  the  higher  current  model 
in  rectifier  substations. 
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The  Transfoshunt' s mode  of  operation  is  shown  in  Figure  9.1.  The 
Transfoshunt  employs  a toroidal  iron  core  with  air  gap.  The  conductor 
carrying  current  to  be  measured  passes  through  the  hole  in  the  toroid,  and 
serves  as  a 1-turn  primary  winding.  A Hall -effect  sensor  is  positioned  in  the 
air  gap.  Hall  voltage  is  sensed  and  amplified  to  drive  current  through  a 
secondary  winding  wound  on  the  toroid,  with  a polarity  that  tends  to  drive  the 
magnetic  field  in  the  air  gap  to  zero.  Thus,  a negative-feedback  system  with 
large  value  of  loop  gain  is  created.  The  effects  of  nonlinearity  of  the  iron 
circuit  are  minimized  by  the  fact  that  the  dc  flux  level  in  the  iron  is 
approximately  zero. 

Current  in  the  secondary  winding  is  sensed  by  an  external  resistor.  DC 
voltage  across  the  sensing  resistor  is  proportional  to  dc  current.  AC  voltage 
across  the  resistor  is  proportional  to  ac  current,  and  may  be  amplified, 
recorded,  or  observed  on  a spectrum  analyzer.  Figure  9.2  shows  a typical 
instrumentation  arrangement  for  observing  and  recording  conductive 
interference. 

One  advantage  of  the  Transfoshunt  is  that  it  is  a dc  instrument,  and  does 
not  suffer  from  the  low-frequency  cutoff  effects  of  standard  current 
transformers.  Ferrite-core  current  transformers  may  be  used  to  sense 
conductive  interference  currents  at  frequencies  above  the  cutoff  frequency  of 
the  current  transformers  and  at  dc  current  levels  low  enough  to  avoid  core 
saturation. 

9.3  Representative  data  from  one  rapid  transit  car 


Figure  9.3  shows  data  recorded  during  a test  of  a particular 
chopper-controlled  rapid  transit  car  at  the  Transportation  Test  Center  (TTC) 
near  Pueblo,  Colorado.  A 10,000-ampere  Transfoshunt  and  associated 
instrumentation  was  installed  in  Thyristor  Substation  No.  2,  and  data  was 
recorded  as  the  rapid  transit  car  accelerated  past  the  substation  third-rail 
feed  point.  Data  shown  in  Fig.  9.3  was  obtained  by  running  the  spectrian 
analyzer  in  the  "Maximum"  mode,  which  captures  the  largest-magnitude  value  of 
each  spectral  line  over  the  duration  of  a run.  Runs  had  durations  of  appx.  30 
seconds. 
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FIGURE  9.3. 


Spectral  plot  of  conductive  interference  from 
a one-car  train. 
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The  substation  is  a three-phase  full -wave  rectifier  having  a fundamental 
ripple  frequency  of  360  Hz.  The  chopper  car  in  question  had  a chopper 
frequency  of  400  Hz.  Substation  and  car  harmonics  are  seen  to  be  clearly 
identifiable  and  measureable  throughout  the  audio-frequency  signaling 
spectrum. 

The  thyristor  substation  employed  at  TTC  is  not  a standard  type  used  in 
rapid  transit  systems.  For  a description  of  the  audio-frequency  harmonic 
characteristics  of  standard  diode  rectifier  substations,  see  Nene  [Ref.  13]. 

9.4  Conductive  Interference  Data  From  Multi-Car  Trains 


Measurements  of  chopper-induced  conductive  interference  currents  flowing 
on  the  third  rail  have  been  made  for  trains  of  1 to  6 cars  in  length.  These 
measurements  were  made  on  the  BART  high-speed  test  track  in  Hayward,  CA. 
Figure  9.4  shows  the  configuration  of  rails,  dc  power  feed  from  substation, 
and  instrumentation  used  in  this  test.  Measurements  were  made  with  trains 
accelerating  from  rest  away  from  the  measurement  point,  and  with  trains 
accelerating  from  rest  away  from  a point  800'  (246  meters)  distant  from  the 
measuring  point. 

A spare  BART  car  was  parked  at  the  measuring  point  to  house  the 
instrumentation  and  provide  power  and  air  conditioning.  The  line  filter 
reactor  of  this  car  was  jumpered  around,  so  that  the  car's  line  filter 
capacitor  bank  served  as  a 12,000  mf  shunt  from  third  rail  to  dc  return  at  the 
measuring  point.  This  procedure  removed  the  inductance  of  the  substation  and 
dc  feeder  cables  from  the  network,  and  created  two  separate  loops  that  met  at 
the  capacitor.  Substation  harmonics  flowed  in  one  loop,  and  chopper  harmonics 
in  the  other.  Two  10,000-ampere  Transfoshunts  were  used  to  separately  observe 
substation  and  chopper  harmonics.  Spectral  plots  similar  to  those  of  Figure 
9.3  were  obtained. 

To  observe  the  behavior  of  chopper  conductive  interference  currents  as  a 
function  of  time,  the  tape-recorded  signals  from  the  third-rail  Transfoshunt 
were  played  back  through  a swept-frequency  spectrum  analyzer  set  up  to  serve 
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FIGURE  9.4.  Circuit  for  measuring  multi-car 
conductive  interference. 
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as  a narrow-band  filter  and  detector.  The  spectrum  analyzer  was  adjusted  to 
zero  Hz  sweep  width,  30  Hz  resolution  bandwidth,  and  5 seconds/division  sweep 
speed.  The  vertical  scale  was  linear,  and  was  calibrated  by  using  a known 
reference  signal  injected  by  means  of  a multi -turn  reference  winding  on  the 
Transfoshunt. 

Data  are  shown  for  the  fundamental  chopper  harmonic  of  218  Hz,  and  the  7th 
chopper  harmonic  at  1526  Hz.  Note  that  the  vertical  scale  differs  from 
picture  to  picture. 

Figure  9.5a  shows  218  Hz  data  for  one  car  starting  from  rest  at  the 
measuring  point.  Note  that  harmonic  current  increases  and  then  decreases  as 
the  chopper  duty  cycle  starts  from  a low  value  and  increases  to  nearly  100%. 
Maximum  harmonic  content  is  reached  when  duty  cycle  is  50%.  Ten  seconds  after 
starting,  there  is  a step  increase  in  current  when  the  field  shunting  occurs. 

Comparing  Figures  9.5a-f  for  1 to  6 cars  respectively,  one  observes  a 
random  behavior  of  peak  amplitude  with  number  of  cars,  as  theory  predicts.  In 
studying  this  data,  one  should  remember  that  over  the  30-second  interval  in 
each  run  for  which  significant  harmonic  current  exists,  the  relative  phases  of 
interference  currents  from  different  cars  do  not  have  a chance  to  change  with 
respect  to  each  other.  However  over  the  minutes  required  prepare  for  the  next 
run,  the  phases  do  have  time  to  re-randomize. 

In  observing  the  effect  at  218  Hz  of  moving  the  starting  point  from  0 to 
800'  away  for  a specific  number  of  cars,  one  observes  the  general  tendency  for 
the  increased  distance  to  have  a far  greater  effect  in  reducing  peak  current 
for  the  longer  trains  than  for  the  shorter.  This  is  once  again  in  keeping 
with  the  analytical  results  outlined  in  Section  7. 

Figures  9.5g  and  h,  taken  at  1526  Hz  for  1 and  2 cars,  demonstrate  clearly 
an  additional  dynamic  effect— the  periodic  variation  in  strength  of  a given 
harmonic  as  a function  of  duty  cycle.  The  effects  of  increasing  the  number  of 
cars,  and  of  moving  the  starting  point,  are  seen  to  be  generally  the  same  at 
1526  Hz  as  at  218  Hz. 
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FIGURE  9.5.  Conductive  interference  current  at 

first  and  seventh  chopper  harmonics 
(218  and  1526  Hz)  for  BART  trains. 
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FIGURE  9.5,  cont'd 
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FIGURE  9.5,  cont'd 
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Table  9.1  summarizes  the  peak  currents  observed  for  a set  of  single  runs 
of  trains  with  1 to  6 cars.  Data  are  shown  for  first  and  seventh  harmonics, 
for  the  two  distances.  Remember  when  reading  this  data  that  all  values  for 
multi -car  trains  must  be  regarded  as  specific  outcomes  of  measurements  of 
random  quantities.  The  data  demonstrate  the  randomness  of  the  process,  but 
also  demonstrate  the  effect  of  added  circuit  inductance  in  decreasing  Cl 
current.  The  fact  that  currents  do  not  increase  dramatically  with  train 
length  is  also  evidenced. 


TABLE  9.1 


Peak  third-rail  current  in  amperes  rms  for  first  and 
seventh  chopper  harmonics  for  multi -car  trains  and  for 
two  circuit  lengths. 


218  Hz  1526  Hz 


cars 

0 m 

244  m 

0 m 

244  m 

1 

37  a 

27  a 

0.20  a 

0.17  a 

2 

60 

10 

0.38 

0.24 

3 

80 

30 

0.38 

0.24 

4 

47 

22 

0.27 

0.15 

5 

35 

6 

0.40 

0.22 

6 

38 

13 

0.38 

0.25 

9.5  Summary 


The  techniques  for  measuring  Cl  currents  in  the  third  rail  or  conductors 
feeding  the  third  rail  are  well  documented  and  have  been  used  on  numerous 
occasions.  Very  repeatable  results  for  chopper  harmonic  current  amplitudes 
can  be  obtained  for  single  rapid  transit  cars.  However,  chopper  harmonic 
current  levels  from  multi -car  trains  are  statistically  distributed,  and  data 
vary  statistically  from  run  to  run. 
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Rigorous  correlation  of  field  data  with  the  statistical  results  of  Sec.  8 
would  require  the  collection  of  much  more  data--enough  to  provide  a 
statistical  data  base.  This  data  would  have  to  be  collected  under  the  same 
carefully  controlled  conditions  as  detailed  in  Sec.  9.4. 

Due  to  the  cost  and  complexity  of  making  measurements,  it  was  only 
possible  to  make-  one  good  run  from  each  starting  location  for  each  train 
length.  However,  the  results  of  the  data  do  fit  in  a general  fashion  with 
what  is  predicted  by  the  circuit  analysis  and  statistics  outlined  in  previous 
sections  of  this  report. 
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10.  CONCLUSION 


Conductive  interference  in  rapid  transit  signaling  is  a problem  that  is 
complex,  due  to  the  physically  large  nature  of  circuits,  and  the  number  of 
separate  circuit  elements  involved.  In  addition,  statistical  considerations 
must  be  taken  into  account.  However,  the  individual  aspects  of  conductive 
interference  are  tractable  as  far  as  analysis  and  testing  is  concerned. 

We  have  examined  different  aspects  of  this  problem  and  the  contributions 
of  many  workers  from  a common  perspective.  Circuit  models  have  been  developed 
for  circuits  comprised  of  track  and  third  rail  using  standard  transmission 
line  and  lumped-circuit  techniques. 

Recommended  practices  have  been  developed  for  observing  and  measuring  Cl 
currents,  and  the  techniques  contained  in  these  have  been  used  to  collect  data 
contained  in  this  report. 

The  efforts  reported  here  were  made  in  response  to  specific  Cl  problems 
that  arose  in  new  rapid  transit  systems.  The  results  that  have  been  obtained 
should  assist  in  the  avoidance  or  solution  of  similar  Cl  problems  in  the 
future.  For  instance.  Section  6 showed  that  in  the  case  of  balanced  two-rail 
track  circuits,  the  injection  of  Cl  currents  into  the  receiver  leads  is  de- 
pendent on  the  asymmetrical  location  of  the  third  rail  with  respect  to  the 
running  rails.  One  way  to  avoid  this  potential  problem  would  be  the  use  of 
overhead  catenaries  instead  of  third  rails.  This  is  but  one  example.  In  the 
design  of  rapid  transit  systems,  every  single  aspect  of  design  should  be 
examined  at  the  outset  for  its  impact  on  Cl. 
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APPENDIX 

FORTRAN 

PROGRAM  FOR  CALCULATING  THE 

CONDUCTIVE 

INTERFERENCE  TRANSFER  FUNCTION 

ZT,  S,  Zl,  Z2,  YPRIME , HCOND 
ZRPARZl 


00100C************************************************************************* 
00110  PROGRAM  TRANS (INPUT, OUTPUT) 

00120C************************************************************************* 

00130C 

00140C 

00150C  *****TYPE  DECLARATIONS***** 

00160C 

00170C 

00180  COMPLEX  ZSER , YSH,  ZO,  GAMMA,  ZR, 

00190  COMPLEX  PHASE,  CSH,  CCH,  ZTPARZl, 

00200  COMPLEX  CXP,  CXM,  CCHM1 , PHSQ 

00210  REAL  M3RR,  LTK 

00220C 

00230  PRINT* , 1 *^*^Tfc,****^*******^*^*^***vt******,**nte**^****^*******'*-'^f •*■■**••*'■*■•*•■*■'*'■*'■*■*■*'  * 
00240  PRINT*,'  ' 

00250  PRINT*,'  *****PROGRAM  TRANS*****' 

00260  PRINT*,  ' ' 

00270  PRINT*,' 

00280  PRINT*,' 

00290  PRINT*,' 

00300  PRINT*,' 

00310  PRINT*,' 

00320  PRINT*, ' 


THIS  PROGRAM  COMPUTES  THE  COMPLEX  TRANSFER  FUNCTION  GIVING' 
THE  RATIO  OF  RECEIVER  INPUT  CURRENT  IN  A TRACK  CIRCUIT  TO' 
THE  CURRENT  IN  THE  THIRD  RAIL.  IT  IS  WRITTEN  IN  FORTRAN  77,' 
FOR  INTERACTIVE  USE.  DATA  INPUT  PROMPTS  ARE  GIVEN  ON  THE' 
SCREEN,  AND  OUTPUT  IS  TO  THE  SCREEN  OR  TO  A PRINTER  DRIVEN' 
OFF  THE  VIDEO  TERMINAL.' 

00330  PRINT*, ' ' 

00340  PRINT*,'  irirkicicicicieicSckie-kicicirieickirkirkickicicirickirick-kieicidi-kidckirkiriiirirkirk-irirki^-k-k-kicirk-k-k-k  * 

00350  PRINT*, ' ' 

00360  50  CONTINUE 

00370C 

00380  PRINT*,'  *****DATA  INPUT*****' 

00390  PRINT*,'  ' 

00400  PRINT*, 'TO  CONTINUE  ENTER  TRACK  LENGTH  D GREATER  THAN  ZERO.' 

00410  PRINT*, 'TO  STOP  ENTER  TRACK  LENGTH  D EQUAL  ZERO.' 

00420  PRINT*, 'ENTER  TRACK  CIRCUIT  LENGTH  D IN  METERS  — D =' 

00430  READ* ,D 
00440C 
00450C 
00460C 
00470  IF 
00480 
00490 
00500 
00510 
00520 
00530 
00540 
00550 
00560 
00570 
00580 
00590 
00600 


*****BEGIN  the  main  if  block***** 

(D.GT.0.01)  THEN 
PRINT* , ' ' 

PRINT*, 'TYPICAL  VALUES  OF  THIRD  RAIL  TO' 
PRINT*, 'RUNNING  RAIL  IMPEDANCE  M3RR  ARE' 


MARTA 

NYCTA 

WMATA 

MBTA 


0.234  UH/METER' 
0.234  UH/METER' 
0.241  UH/METER' 
0.273  UH/METER' 


PRINT*, ' 

PRINT*, ' 

PRINT*, ' 

PRINT* , ' 

PRINT* , ' ' 

PRINT*, 'ENTER  M3RR  IN  UH/METER  — M3RR  = ' 

READ* ,M3RR 

PRINT*, 'ENTER  SERIES  TRACK  RESISTANCE  RTK  IN  UOHMS /METER  — RTK  = ' 
READ* , RTK 

PRINT*, 'ENTER  SERIES  TRACK  INDUCTANCE  LTK  IN  UH/METER  — LTK  = ' 


-A.l  - 


00610 

00620 

00630 

00640 

00650 

00660 

00670 

00680 

00690 

00700 

00710 

00720 

00730 

00740 

00750 

00760 

00770 

00780 

00790C 

00800C 

00810C 

00820C 

00830C 

00840 

00850 

00860 

00870 

00880 

00890 

00900 

00910 

00920 

00930 

00940C 

00950 

00960 

00970 

00980 

00990 

01000 

01010 

01020 

01030 

01040 

01050C 

01060 

01070 

01080C 

01090 

01100 

oino 

01120 

01130 

01140C 


READ*,  LTK 

PRINT*, 'ENTER  REAL  AND  IMAGINARY  PARTS  OF  ZT  — ZTR,ZTI  = ' 

READ*,  ZTR,  ZTI 

PRINT*, 'ENTER  REAL  AND  IMAGINARY  PARTS  OF  ZR  — ZRR,ZRI  = ' 

READ*,  ZRR , ZRI 

PRINT*, 'ENTER  BALLAST  IN  OHM-KFT  — BALLAST  = ' 

READ*,  BALLAST 

PRINT*, 'ENTER  FREQUENCY  IN  HZ  — F = ' 

READ*,  F 
PRINT* , ' ' 

PRINT*, 'TYPE  PARAMETER  K = 0 FOR  BALANCED  2-RAIL  TRACK  CIRCUIT' 
PRINT*,-'  = 1 FOR  TK  CKT  WITH  SIGNAL  RAIL  NEXT' 

PRINT* , ' TO  THIRD  RAIL ' 

PRINT*,'  = -1  FOR  TK  CKT  WITH  DC  RETURN  RAIL  NEXT' 

PRINT*,'  TO  THIRD  RAIL' 

PRINT*,'  ' 

PRINT*, 'ENTER  TYPE  PARAMETER  K — K=  ' 

READ*,  K 


*****C(3MP  UTAT I ON***** 


YSHR  = 1 . / (304.7*BALLAST) 

YSHI  = 0. 

YSH  = CMPLX (YSHR, YSHI) 

ZSER  = CMPLX ( RTK* 1 .E-6,  6 .28319*F*LTK*1 .E-6 ) 

ZT  = CMPLX (ZTR, ZTI) 

ZR  = CMPLX (ZRR, ZRI) 

ZO  = CSQRT( ZSER /YSH) 

GAMMA  = CSQRT(ZSER*YSH) 

PHASE  = GAMMA*D 

PHMAG  = SQRT (REAL (PHASE )**2  + AIMAG(PHASE)**2) 

IF  (PHMAG .LE.0.01)  THEN 
PHSQ  = PHASE**2 
CSH  = PHASE* ( 1 . + PHSQ/6.) 

CCHMl  = (PHSQ/ 2 . )*( 1 . + PHSQ/12.) 

ELSE 

CXP  = CEXP (PHASE) 

CXM  = CEXP (-PHASE) 

CSH  = (CXP  - CXM)  / 2 . 

CCHMl  = ((CXP  + CXM)/2.)  - 1. 

END  IF 

Z1  = Z0*CSH/CCHM1 
Z2  = Z0*CSH 

IF  ((ZT.EQ.O.).OR.(Zl.EQ.O.))  THEN 
ZTPARZ1  = 0. 

ELSE 

ZTPARZl  = ZT*Z1/(ZT  + Zl) 

END  IF 
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01150 

01160 

01170 

01180 

01190 

01200C 

01210 

01220 

01230 

01240 

01250 

01260 

01270 

01280 

01290C 

01300C 

0131 0C 

01320C 

01330C 

01340 

01350 

01360 

01370 

01380 

01390 

01400 

01410 

01420 

01430 

01440 

01450 

01460 

01470 

01480 

01490 

01500 

01510 

01520 

01530 

01540 

01550 

01560 

01570 

01580 

01590 

01600 

01610 

01620 

01630 

01640 

01650 

01660 

01670 

01680 


IF  ((ZR.EQ.O) .OR.(Zl.EQ.O))  THEN 
ZRPARZ1  = 0. 

ELSE 

ZRPARZl  = ZR*Z1/(ZR  + Zl) 

END  IF 

YPRIME  = (Z1/(Z1  + ZR) ) / (ZTPARZl  + Z2  + ZRPARZl) 

S = CMPLX(0 . , 2.*3.14159*F) 

AK  = REAL(K) 

HCOND  = YPRIME* (S*M3RR*1 .E-6  + (AK*ZSER/2 . ) )*CSH/GAMMA 
HCONDR  = REAL (HCOND) 

HCOND I = AIMAG(HCOND) 

HCONDM  = SQRT (HCONDR**2  + HCONDI**2) 

HCONDPH  = ATAN(HCONDI/HCONDR)*180./3. 14159 


***** DATA  OUTPUT***** 


= ' ,G11 .4, ' UH/M') 

= ' ,G11 .4, ' + J( ' ,G11 .4, ' ) OHMS/M') 
= ' ,G11 .4, 1 UH/M') 

= ' ,G11 .4, ' + J( ' ,G11 .4, ' ) OHMS') 

= ' ,G11 .4, ' + J( ' ,G11 .4, ' ) OHMS') 


PRINT*  , 1 ■*"*'5ff'*,,*"*-'*,,*"*,,A"5Sr-*"*,'*,'5Sc,'*"*,,5fe'‘*’-5fe,,*“*,,*-'*,,*'*********************************'*,*r  • 

PRINT*,'  ' 

PRINT*,'  *****iNPUT  DATA*****' 

PRINT*,'  ' 

PRINT  11,  D,  F,  K 

11  F0RMAT(1X,'D  = *,G11.4,'  M,  F = ',G11.4,'  HZ,  K = ',12) 
PRINT*,'  ' 

PRINT  12,  M3RR 

12  FORMAT ( IX, 'M3RR 
PRINT  13,  ZSER 

13  FORMAT (IX, 'ZSER 
PRINT  14,  LTK 

14  FORMAT ( IX, 'LTK 
PRINT  15,  ZT 

15  FORMATQX,  'ZT 
PRINT  16,  ZR 

16  FORMAT ( IX, 'ZR 
PRINT  17,  BALLAST,  YSHR 

17  FORMAT (IX, 'BALLAST  =' ,G11 .4, 'KHOM-FT , YSHR  =' ,G11 .4, 'MHOS/M' ) 
PRINT* , ' ***********************************************************  1 
PRINT*,'  ' 

PRINT* , ' *****0UTPUT  DATA***** ' 

PRINT*,'  ' 

PRINT  18,  PHASE 

18  FORMAT ( IX , ' GAMMA*D  = ',G11.4,'  + J( ' ,G11.4, ' ) ' ) 

PRINT  19,  ZO 

19  F0RMAT(1X, 'ZO  = ',G11.4,'  + J(',G11.4,')  OHMS') 

PRINT  21,  Zl 

21  F0RMAT(1X, 'Zl 
PRINT  22,  Z2 

22  FORMAT (IX , 'Z2 
PRINT  23,  HCOND 

23  FORMAT (IX, 'HCOND 
PRINT  24,  HCONDM,  HCONDPH 

24  FORMAT (IX, ' 


= ' ,G11 .4, ' + J( ' ,G11 .4, ' ) OHMS') 

= ' ,G11 .4, ' + J( ' ,G11 .4, ' ) OHMS') 

= ' ,G11 .4, ' + J( ' ,G11 .4, ' ) ' ) 

= ' ,G11 .4, ' AT  ' ,G11 .4, ' DEGREES') 
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01690  PRINT* , ' *,***'***********,**'************'***'**,*********'***'*,******,****ifr*'*  • 

01700  PRINT*,'  ' 

01710  PRINT*,'  ' 

01720C 

01730C  *****PROGRAM  WILL  NOW  LOOP  TO  REPEAT  FOR  NEW  DATA***** 

01740  GO  TO  50 

01750C 

01760  END  IF 
01770C 

017 80C  *****THIS  ENDS  THE  MAIN  IF  BLOCK***** 

01790C 

01800  PRINT*  , ' *******'*'**'*'***'***'****'*'*,*,*'***,**********'****'***'******'*■*'*** 1 
01810  PRINT*, ' ' 

01820  PRINT*,'  COMPUTATION  SESSION  CLOSED.  ADIOS,  AMIGO.' 

01830  PRINT*, ' ' 

01840  PRINT* , ' ************■***'*'****************'****,*'***,***'*'********,**■*•* ' 

01850  END 
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